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SUMMARY 


Snow-drifting  characterlatic*  of  various  structurea  were  studied  by 
scale-model  wind-tunnel  experiments.  Modeling  criteria  developed  by 
two  different  methods  were  used  as  a  basis  for  selecting  a  material  for 
simulating  snow  in  scale-model  experiments  and  for  determination  of 
test  variables. 

A  number  of  materials  were  tested  for  physical  properties  which 
would  make  them  suitable  for  snow  simulators.  It  was  found  that,  for  a 
1/10  model  scale,  0.  01  cm  diam  crystalline  borax  was  the  beet  available 
material. 

The  work  performed  in  the  3^  by  7  ft  wind  tunnel  adequately  demon¬ 
strated  the  feasibility  of  using  geometrically  and  physically  scaled  syn¬ 
thetic  snow  in  a  wind  tunnel.  It  was  also  shown  that  it  is  possible  to 
present  several  years  of  arctic  snowstorm  environment  in  a  few  hours 
of  wind-tunnel  time.  The  actual  tests  showed  these  major  points: 

1 .  Close  spacing  of  buildings  will  result  in  coalescence  of 
drifts.  This  applies  regardless  of  the  orientation  of  the  long  axis  of  the 
building  with  respect  to  wind  flow. 

2.  When  rectangular  buildings  must  be  grouped  in  closely 
spaced  units,  it  appears  that  they  should  be  erected  with  the  long  axis 
parallel  to  the  dominant  wind  direction.  Channeling  the  wind  flow  be¬ 
tween  paired  buildings  oriented  with  the  wind  tends  to  keep  the  space 
between  the  longitudinal  rows  of  buildings  open. 

3.  Snow  erosion  may  occur  beneath  buildings  erected  on  columns 
above  the  surface.  The  columns  must  be  carried  to  an  unerodable  stratum 
or  they  must  be  set  on  a  crib  foundation  to  prevent  undermining  of  the 
footing. 

4.  A  V-shaped  snow  fence  produced  a  clear  area  downwind  for 
a  distance  equal  to  approximately  25  times  the  height  of  the  fence.  This 
phenomenon  may  be  applied  to  the  protection  of  trench  entrances  to  under¬ 
snow  structures. 

5.  The  tests  with  the  borax  snow  simulator  produced  classic 
patterns  of  drifting  very  similar  to  those  characteristic  of  cold  dry  arctic 
snow. 

6.  The  drift  patterns  produced  by  the  unidirectional  flow  in  the 
tunnel  probably  were  more  concentrated  longitudinally  than  those  developed 
under  natural  conditions  where  some  variation  in  wind  direction  and  speed 
occurs  during  and  following  periods  of  blowing  and  drifting  snow. 
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SCALE  MODEL  STUDIES  ON  SNOW  DRIFTING 

by 

Gordon  H.  Strom,  George  R.  Kelly, 

Edwin  L.  Keitz  and  Robert  F.  Weiaa 

INTRODUCTION 

One  of  the  more  troubleaome  featurea  of  arctic  life  haa  been  the  aerioua  handicap 
impoaed  on  all  forma  of  activity  by  blowing  and  drifting  anow.  Over  the  yeara  man  haa 
developed  techniquea  to  gain  partial  control  of  thia  phenomenon.  Such  devicea  aa  anow 
fencea  probably  were  evolved  from  obaervation  of  the  effect  of  auch  natural  windbreaka 
aa  hedgerowa  on  anow  drifting  and  owe  their  deaign  more  to  local  availability  of  materi- 
ala  than  to  aerodynamical  or  meteorological  conaiderationa.  The  anow  fencea  uaed  in 
varioua  countriea  differ  conaiderably  in  form.  It  ia  difficult  to  account  for  the  particular 
configuration  which  ia  in  local  uae  except  that  it  ia  traditional  and  low  in  coat. 

One  of  the  reaaona  for  the  abaence  of  more  conacioua  deaign  in  the  development  of 
anowdrift  control  methoda  haa  been  the  neceaaity  and  difficulty  of  field  teating  any  newly 
conceived  model  over  a  period  of  yeara  in  order  to  obaerve  ita  performance  over  a  full 
range  of  meteorological  conditiona.  However,  a  few  decadea  ago,  the  idea  of  uaing  a 
wind  tunnel  to  teat  auch  deaigna,  with  conaequent  aavinga  in  time  and  coat,  occurred  to 
people  in  aeveral  anow  countriea.  E.  A.  Finney  (1940)  in  Michigan  performed  a  variety 
of  wind  tunnel  experimenta  uaing  balaa  aawduat  and  mica  flakea  to  repreaent  anow. 

From  theae  he  drew  a  number  of  intereating  concluaiona  about  the  aizea  of  turbulent 
regiona  in  relation  to  the  obataclea  which  produced  them.  The  reaulta  of  hia  atudiea 
have  been  found  very  uaeful  by  people  working  on  highway  deaign  and  maintenance. 

In  Scandinavian  countriea  alao  there  waa  intereat  in  wind-tunnel  experimentation  on 
blowing  anow.  C.  Nokkentved  (1938)  performed  aimilar  experimenta  and  obtained  re- 
aulta  which  were  in  general  accord  with  Finney'a  findinga.  Becker  (1944)  conducted  a 
few  teata  in  Germany,  but  later  deapaired  of  obtaining  entirely  aatiafactory  reaulta 
from  auch  inveatigationa  in  view  of  the  aeemingly  inaurmountable  problema  of 

(a)  aimulating  the  anow  adequately,  and 

(b)  aatiafying  the  acaling  criteria  accurately. 

S.  Hallberg  and  othera  (1943)  of  the  State  Road  Inatitute  of  Sweden  inveatigated  the 
eddiea  around  anow  fencea,  uaing  full-acale  fencea  in  the  out-of-doora  with  ailk 
atreamera  to  define  the  flow  patterna. 

Theae  inveatigationa  proved  to  be  quite  uaeful  in  eatabliahing  certain  groaa  concepta 
of  anow  blowing  and  drifting  behavior,  but  all  of  them  auffered  from  the  difficultiea 
noted  by  Becker.  In  general,  the  materiala  uaed  to  aimulate  anow  were  choaen  arbi¬ 
trarily,  baaed  upon  auch  aecondary  conaiderationa  aa  coat,  availability,  or  auperficial 
likeneaa  to  actual  anow.  Although  the  neceaaity  for  acaling  the  atructurea  repreaented 
waa  recognized,  little  waa  done  to  inaure  that  Ae  particlea  of  blowing  anow  were  acaled 
in  proportion.  Conaequently,  the  end  reault  waa  a  little  like  examining  the  behavior  of 
flying  baaeballa  aa  a  baaia  for  drawing  concluaiona  on  the  behavior  of  blowing  anow. 

Or.  R.  W.  Gerdel  (1960]tof  U.  S.  Army  Snow,  Ice  and  Permafroat  Reaearch 
Eatabliahment ,  *  pointed  out  that  it  ia  eaaential  to  acale  the  dimenaiona  and  phyaical 
propertiea  of  the  material  aimulating  anow  in  wind  tunnel  inveatigationa  in  order  to  gain 
true  aimilitude  to  the  behavior  of  the  anow. 

The  preaent  project  reaulted  from  thia  auggeation  and  had  aa  ita  goal  the  inveati- 
gation  of  the  micro-meteorological  factora  which  control  the  drifting  of  anow,  with  the 
apecial  objective  of  determining  a  material  which  can  be  uaed  in  the  wind  tunnel  to 
repreaent  the  behavior  of  anow  accurately. 

71[e3eaIgn«te9TJ7~37~Xrmy^T!roI3~l^gIoni~K»earch  and  Engineering  Laboratory, 

1  February  1961. 
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SCALE  MODEL  STUDIES  ON  SNOW  DRIFTING 


THEORETICAL  ANALYSIS  OF  SALTATION  PHENOMENA  IN  SNOW  DRIFTING 

The  details  of  snowdrifting  phenomena  are  so  complex  that  it  is  difficult  to  obtain 
an  understanding  of  the  fundamental  mechanism  involved  by  either  theoretical  or  experi¬ 
mental  methods.  A  combination  of  the  theoretical  and  experimental  approaches  is 
likely  to  yield  more  progress  than  either  one  alone.  The  following  section  presents  a 
theoretical  analysis  of  an  important  phase  of  snow  drifting  and  adds  to  the  understanding 
of  the  fundamental  mechanism. 

Pioneering  work  in  drifting  problems  has  been  carried  out  by  Bagnold  (1941)t  who 
provides  basic  descriptions  of  the  various  significant  phenomena.  The  transport  of 
small  particles  over  a  surface  has  been  shown  to  consist  of  suspension,  saltation,  and 
surface  "creep"  of  the  material  involved.  In  the  application  of  drifting  problems  to  snow, 
particle -surface  interactions  are  of  primary  interest.  The  larger  part  of  particle  trans¬ 
port  over  a  surface  has  been  shown  to  be  due  to  saltation,  or  bouncing  of  particles  from 
point  to  point  in  the  general  direction  of  the  free  stream.  Hence,  it  is  logical  to  first 
consider  the  saltation  process  and  attempt  to  obtain  reasonable  results  for  a  specific 
physical  model.  Such  results  should  include  a  predicted  "wave  length"  of  drifting  and 
vertical  particle  distribution  in  the  saltation  layer. 

The  analysis  consists  of  three  fundamental  subjects:  (1)  analysis  of  the  movement 
of  a  particle  through  the  air,  between  interactions  with  the  surface,  (2)  description  of  a 
particle  field  on  the  surface,  and  interaction  of  an  airborne  particle  with  the  field;  and 
(3)  a  correlation  between  these  phenomena  that  will  provide  a  steady-state  solution  to  the'' 
particle  distribution. 

Development  of  theory 

Assumptions. 

1.  Particles  are  elastic,  spherical,  equally  sized,  and  much  smaller  than  the 
boundary  layer  thickness. 

2.  Field  consists  of  a  special  closely  packed  geometrical  arrangement  of  particles, 
free  to  move  in  the  plane  of  the  field  only. 

3.  Phenomena  are  investigated  at  points  sufficiently  far  downstream  to  allow  a 
large  number  of  particles  to  reach  "equilibrium"  with  the  field. 

4.  Only  primary  interactions  are  considered,  secondary  "ejections"  being  shown 
to  comprise  a  small  fraction  of  the  total  saltation. 

5.  Rolling  particles,  randomly  located,  initiate  saltation  at  all  points  on  the  field. 

6.  In-air  collisions  of  particles,  boundary  layer  turbulence,  and  crosswind  are 
neglected. 

Trajectory  analysis.  Figure  1  shows  the  relationship  between  velocities  and  forces 
on  a  particle  at  any  time  during  its  trajectory  and  Figure  2  shows  the  simplified  boundary 
layer  employed  to  facilitate  solution  of  equations  of  motion.  V  is  a  linear  function  of  ^ 
within  the  boundary  layer  depth  The  assumption  is  valid  IT  the  particle  does  not  pass 
above  and  is  generally  indicative  of  a  laminar  sublayer  to  the  turbulent  atmospheric 
boundary  layer. 

The  following  two  equations  are  formed  by  setting  the  sums  of  forces  equal  to  parti-  . 
cle  mass  times  its  accelerations  in  the  x  and  ^  directions.  ^  is  taken  vertically  upward. 
Dg  is  the  aerodynamic  drag  force  which  acts  in  the  direction  opposite  to  relative  velocity 
^  of  the  particle  and  surrounding  air.  mg  is  particle  weight. 


Dg  sin  ^  -  mg  =my 

(1) 

Dg  cos  p  =  mit. 

(2) 

The  velocities  are  related  to  the  angle  ^  by  the  following  equation  where  ^  is  the  vertical 
component  of  particle  velocity: 

sinp=  -y/Vr. 
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Figure  1.  Forces  and  velocities  Figure  2.  Linear  velocity  profile 

of  a  particle.  assumed  for  boundary  layer. 

At  the  low  Reynolds  numbers  encountered  in  these  particle  motions,  the  drag  force  is 
approximately  proportional  to  relative  velocity  as  follows; 

Dg  =  '‘sV,.  (3) 

The  left  side  of  eq  1  becomes 

D^  sinp  =  -koV^y/V^  =  -koir 

and  eq  1  becomes 

y  +  (ko  /m)y  +  g  =  0.  (4) 

In  eq  4,  taking  the  Laplace  transform  of  both  sides, 

s*f(s)  -y(0)  +  ^s  f(s)  +  J  =  0 
which  is  simplified  to: 


taking  the  inverse  transform: 


(5) 
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Equation  2  is  no>v  written 
cos  fi  =  mS 

where 

V-  it 

cos  S  =  ; 

r 

therefore, 


=  mil 


or 

it  +  (ko/m)x  =  (ko/m)£,(t) 
since  V_  can  now  be  expressed  as  f|(t). 

In  operator  form,  eq  6  becomes: 

|d*  +  (ko/m)D  jx  =  (ko/m)f|(t) 
or 

d|d  +  (ko/m)jx  =  (ko/m)f|(t). 

Letting  |d  +  (ko/m)jx  =  H;  DH  =  (ko/m)fi(t). 

Substituting 

k,y  =  V(y)  =  fi(t). 


DH  =  k,[y(0)  +  mg/ko  jj^l  -  -gk,t 

which,  when  integrated,  yields 

[^y(O)  +  mg/ko  j  j*  [*  ■ 


H  =  k, 


^-tko/m)t 


and 

(D  +  ko/m)x  =  k|  jVco)  +  mg/koj|t  =  (m/ko)e  -gk|t*/2  +  C|. 

This  is  a  linear  equation,  having  the  solution: 

X  =  J  e<'''>/*">‘{k,[y(0)  +  mg/k,][t  +  (m/ko)e-<’^« + 

+C,  -k|gt»/2  j-  dt  +  C, 

Integrating  eq  8a  by  parts  gives: 

*  =  ‘‘lA  +  0-  gij  "-i* [Bf  +  -5  ■  llj  +  C,  e  +  -^ 

where  B  =  ko/m  and  A  =  fy(o)  ♦  il 


(6) 


(7) 


(8a) 


(8b) 


I 


( 


V 
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From  the  boundary  condltiona 

X  =  0,  k  =  ic(0)  at  t  =  0;  we  get,  finally: 


X  =  k|A 


X  2  ^  t*  2tl 
+  +  1-  1*^ 


m 


which,  when  combined  with  eq  5,  gives  a  parametric  description  of  the  position  of  a 
particle. 

Leec 

and  form  the  derivatives  of  ^  and  x  with  ^  as  follows 


dx  k|A 
e[F  ^  B 


e'®‘(l  -Bt)  +1  -  ^  -BCe‘ 


Bt 


.-Bt  ,£ 


a?-Ae-  B 

and  the  angle  of  inclination  of  the  particle  velocity,  )|<>  is  given  by: 

•Bt 


Ae 


>g/B 


¥[< 


e'®‘  (1  -  Bt)  +  1 


(10) 


■BCe 


-Bt 


This  quantity  will  be  useful  in  the  following  interaction  analysis. 


Particle  —  field  interaction.  The  most  important  assumption  here  concerns  the 
"model"  snowfield  to  be  used.  For  simplicity  of  analysis  (with  good  approximation  to 
the  physical  situation),  a  "closely  packed"  field  is  assumed  (Fig.  3).  Such  a  field 
seems  quite  arbitrarily  oriented  as  far  as  an  approaching  particle  is  concerned,  and  it 
is  also  reasonable  to  assume  that  field  particles,  though  close  to  each  other,  do  not 
touch  in  general.  The  field  particles  are  then  free  to  move  in  the  plane  of  the  surface, 
but  are  assumed  to  meet  infinite  resistance  normal  to  the  surface.  Furthermore,  an 
approaching  particle  has  two  important  restrictions  on  its  interaction  with  the  field: 

a.  The  maximum  angle  of  approach  off  a  stationary  particle's  centerline  (in  trans¬ 
verse  plane)  is  limited  to  30*  (Fig.  4).  Particles  approaching  as  far  off  center  as 
shown  in  the  figure  will  be  only  slightly  deflected,  because  the  large  transverse  deflec¬ 
tion  angle  is  opposed  by  a  small  transverse  velocity  (compared  to  longitudinal  velocity). 
Thus,  an  essentially  two-dimensional  interaction  may  be  assumed:  i.  e.  ,  spheres 
impacting  on  cylinders  placed  normal  to  the  free  stream. 

b.  The  minimum  distance  between  centers  of  approaching  and  stationary  particles 
is  a  function  of  the  angle  of  approach  (Fig.  5).  This  may  be  expressed  as. 


si.ni(i 


n 

<r 


and  b 


min 


b  . 
min 


o-  (1  -  sintji) 


<r  -  n 


a  -  <r  sinij) 


(11) 


a  relationship  to  be  used  in  actual  calculation. 

Finally,  it  should  be  pointed  out  that  a  typical  occurrence  in  saltation  is  the  ejec¬ 
tion  of  new  particles  from  the  surface.  This  is  indeed  possible  in  the  model  field  set  up, 
as  a  struck  particle  is  free  to  move  until  striking  an  adjacent  particle,  whereupon  it 
will  rebound  into  the  airstream.  Many  more  particles  will,  in  turn,  be  struck,  but  the 


i 
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IV 

TOP  VIEW 


Figure  4.  Maximum  transverse  angle 
of  approach  is  shown  to  be  30*.  (Front 
view  looking  downstream. ) 


suarace 

LSVCR 
SUB  LAYER 

FRONT  VIEW 

(LOOKING  DOWNSTREAM) 

Figure  3.  Assumed  form 
of  particle  field  at  surface. 


Figure  5.  Interaction  limitations. 


loss  of  energy  in  an  interaction  is  so  great  that  even  the  secondary  interaction  results 
in  very  small  ejection  velocities.  In  the  analysis  to  follow,  only  primary  interactions 
are  considered,  but  these  procedures  could  be  extended  to  include  secondary  effects. 

In  the  following  analysis,  equations  are  developed  which  relate  the  velocities  and 
path  angles  of  a  particle  before  and  after  impact  for  one  degree  oi  freedom.  The  im¬ 
perfectly  elastic  properties  are  expressed  in  terms  of  the  coefficient  of  resUtution.  The 
designations  of  velocities  before  and  ^ter  impact  are  shown  in  Figure  6. 

Since  the  total  momentum  of  the  ^proaching  particle's  normal  velocity  component 
may  be  divided  into  horisontal  and  vertical  components,  we  may  say  that  the  change  in 
momentum  of  the  impacting  particle  equals  the  sum  of  the  changes  of  its  components.  In 
this  analysis,  we  of  course  assume  a  smooth  surface  and  smooth  particle.  Hence,  vii 
is  preserved  and  equals  V|  *•  The  vector  diagrams  before  and  after  interaction  are  shown 
in  Figure  6  where  primes  denote  post-collision  quantttiea,  and  p.  is  scattering  angle. 

The  velocities  before  and  after  collision  are  related  to  the  coefficient  of  restitution 
e,  which  is  defined  by: 


The  respective  particles  are  indicated  by  the  subscripts  a  and  b. 

When  applied  to  the  horisontal  components,  the  following  is  the  result: 


KFOUE 


AFTER 


Figure  6.  Vector  diagrame  of  particle  velocity  and 
velocity  components  before  and  after  collision. 


AS  e  -  0.  v,;^  -  v,„^/2.  and  as  e  -  1.  v,'^^  =  0. 

In  the  vertical  direction,  my  is  assumed  infinite  because  of  the  resistance  to  motion 
normal  to  the  surface,  so  that  v*  =  v/  s  0  and 
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When  applied  to  the  vertical  component,  (v^*vg)  and  (vb-va)  are  the  separating  and 
closing  velocities,  respectively,  of  the  two  particles.  Since  there  is  no  resistance  to 
particle  "2"  movement  at  impact,  the  conservation  of  momentum  gives  in  general  form 

V|  +  V,  =  V,"  +  V/  =  VjJ 

v,  =  0  since  particle  "2"  is  initially  at  rest.  The  coefficient  of  restitution  gives 


Solving  for  velocity  V|  in  terms  gives: 


^*nv  "  *'^*nv‘ 

From  Figure  6,  we  may  write: 

''‘nv  ^  ''‘t  ““[j-  -(“+4')] 


tanP.  = 


•  vi'^  cos  j-(a  +  »|»)  +v,-^j^ 

U  •  *" 


where  a  =  sin'*  — ,  as  seen  in  the  figure, 
<r 


and,  of  course,  sin  *(a-t-4<) 


cos 


IT 

2 


■(o  +  lp) 


=  cos  (a  +  il*) 
=  sin  (a  4- 4*)- 


Now,  since  V|  =  v,  cos  a  sin  (a -t-t)*) 

•nv  • 

V|^h  -  '^1  cos  a  cos(o  +  4<) 

and  V|^  s  V|  sina 


(14) 


(15) 


(16) 


and  defining  S 
we  obtain 

tanP^ 


eS  tan  (a  >|i)  *■  S  tan  a 
S  tan  a  tan  (a  4  il*)  1 


(17) 


Particle  distribution.  With  a  uniform  density  of  Iq  particles/ sec -unit  area  im- 
pacting  on  a  field  particle,  the  total  number/sec  passing  through  an  annular  ring  at  a 
distance  b  (see  Fig.  7),  ^  wide  is  dN;  and, 

dN  =  loyb  db 

where  y  =  angular  range  about  center.  (18) 

Since  particles  impacting  at  a  distance  b  rebound  at  the  angle  Pg,  W  is  the  number  of 
particles/sec  rebounding  at  angles  between  Pg  and  Pg  dPg.  We  must  now  compute  b, 
or  a,  and  A  in  terms  of  Pg  and  the  other  relevant  parameters. 

Differentiating  both  sides  of  eq  17  with  respect  to  a,  and  noting  that  d(a  ■)-  4<)  =  we 

get 
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lec*  P  dp  [S  tana  tan(a  +  i|»)  -1  )• 
db  =  — - - = - - - . 

S  •|•ec*  (a  i-  i|<)  [  e  -I-  S  tan*a]  +  sec*a  [  1  +  eS  tan*  (a  +  ( 

Using  tan(a  +  4<)  =  i*"tan a^*anV *  rewrite  eq  17  as: 

tan*a  [S  tan  p  -S  tantp]  -t-  tana  [(S-1)  tanij«  tanp  -t-  S  (e  +  1)]  +  [tan  p  +  eStani|j  =  0 

•  8  8^ 

a  quadratic  in  tana,  which  may  be  solved  to  give  the  required  relationship  for  b  -  b 
(Pg,  etc. ),  i.  e.  , 


tan  a  = 


-B  *  Vb*  -4AC 

- Za - 


and  tana  = 


We  may  now  express  the  number  of  particles  rebounding  per  sec  per  unit  angle  pg  as 

=  K2F  (Pg,  »|/,  e,  <r)  and  call  this  quantity  the  probability  that  a  particle  will  rebound 

at  an  angle  P,,  having  approached  the  field  at  angle  ij<.  Incorporating _K2,  a  constant, 
into  this  probability  (since ^0  is  unknown  anyway),  we  have  obtained 

P  (Ps.  4')  =  f(Pe.  4'.  e,  <r).  (21) 

For  purposes  of  comparing  this  probability  function  with  those  for  the  well  known  cases 
of  infinite  restraint  and  zero  restraint  (in  all  directions),  we  let  4i  =  0  and  e  =  1.  The 
results  are  shown  in  Figure  8. 

Having  obtained  the  trajectory  and  surface  interaction  relationships,  we  attempt  to 
correlate  this  information  to  obtain  the  fundamental  unknown  in  the  problem,  P  (Pg  , 

Vo)i  the  probability  that  a  particle  will  leave  the  surface  at  an  angle  Pg  and  velocity 
Vg.  The  particle  distribution  above  the  surface  may  then  be  calculateo  and  the  "wave 
length"  of  the  saltation  phenomenon  estimated. 

Correlation.  The  following  notation  will  be  used: 

Psg'  ''0  angle  and  velocity  at  beginning  of  a  trajectory  (initial  values) 

Ip,  V|  are  angle  and  velocity  at  end  of  trajectory 

Ps|>  V  angle  and  velocity  after  interaction,  hence  are  initial  values 
of  next  trajectory. 

The  main  assumption  employed  here  is  as  follows:  If  there  is  any,  set  of  probabili¬ 
ties  P(Pgg,  Vq)  that  exist  at  various  values  of  x,  and  are  invariant  inUme,  because  of  a 
boundary  condition  at  Xg ,  we  see  that  the  same  set  exists  at  point  x  +  Ax  because  of  the 
same  boundary  conditions  at  Xg  +  Ax.  The  boundary  conditions  are.  in  fact,  duplicated 
at  X  g  +  Ax  (£or  all  values  of  Ax)  in  a  field  of  randomly  originating  rolling  particles. 

possible  member  of  the  P-  set  are  superimposed  at  all  x's  sufficiently  far 
do^stream.  and  we  can  assume  the  sum  of  these  to  equal  the  P(pr.  vX).  that  is  invari¬ 
ant  with  respect  to  both  time  and  x.  This  argument  reduces  to  the  assumption  that 

P(Pg,.  y{)  =  P(P,^.  vj)  =  P(p^.  V').  (22) 

We  now  look  at  a  set  of  trajectories  terminating  at  a  fixed  point  x.  each  trajectory  being 
determined  uniquely  by  specifying  ip  and  We  then  assume  the“ixistence  of  a  I^iP,  v.f. 
Ttf  *  particle  will  terminate  its  trajectory  at  angle  ip  and  velocity  * 

Then,  since  is  related  to  v,'  and  Pg,.  there  is  a  determinable  number  of  particlS/aec 
^rminating  at  ip  and  a  such  that  the  percentage  rebounding  at  Pg  will  have  velocity  vf. 
The  percentage  rebounding  at  pg.  is.  of  course.  P(Pg, .  iP).  as  already  obtained.  If  wi"^ 
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sum  over-all  values  of  i|i  ,  we  obtain  Uie  total  number/eec  rebounding  at  Pg|  and  W*, 
or  P(Paii  y{).  Finally,  we  note  that  the  number  of  particles  at  t|«  and  V|  equals  num¬ 
ber  having  left  the  surface  at  values  of  ^ae  Xo*  determined  by  relationships  derived 
in  Figure  6.  Thus.  ~ 

pi|i  max 

P(Pa  .Vi")  =  \  P(P  .  vS)  P(p  .  i)<)  d|i  (23) 

‘^ipmin  •«  ■» 

where  6  and  vi  are  functions  of  ii«,  and  vt.. 

•o  •» 

The  required  relationship  between  Vf  and  v{,  mentioned  above,  is  now  derived. 

From  Figure  8  we  see  that: 

(V|')*  =  [v|'^^  +  vi'j  cos  (a  +  +  Ivi'j  sln(i  +  »|»)  + 

and,  with  substitutions  as  before, 

V,'  =  Ei  vi  (24a) 

where  E  =  [ecosu  sin(a  i|«)  -f  sina  cos(a  +  4*)]*  (24b) 

■f  [sin  a  sin(a  +  tf)  cos  a  cos(a 

Finally,  since  V|^  =  V|  sin  4*, 


h'«inPso+-|J  -i 

v»  - - -  (25) 

where  ta  is  the  terminal  time. 

Proar^  for  numerical  solution.  The  program  that  must  be  followed  for  numerical 
s  olution  oi  the  problem  is  now  outlined: 

(1)  and  v/  are  chosen  arbitrarily. 

(2)  a  is  determined  from  eq  20  for  various  41’s. 

(3)  ^  and  ^  are  calculated  from  eq  24a  and  24b  for  various  41' s. 

(4)  y  (0)  is  calculated  from  eq  25  for  various  4''s. 

(5)  d  (0)  is  calculated  from  eq  10  for  various  tj^'s. 

(6)  then  p^  ,  ^  are  easily  calculated  for  various  vj^'s. 

(7)  P(P,  >  4*)  i*  determined  from  eq  21  for  various  4''s. 

(8)  P(P.  ,  Vg)  is  assumed  to  be  of  the  simple  gaussian-type  form: 

•0 


where^  and  C4  are  constants  to  be  determined,  equal  to  the  one-half  the  most  probable 
angle  and*one-half  the  most  probably  velocity  v^g,  respectively;  i.  e. , 

•P  0 


when 


8P 

•0 


=  0, 


8P 

■5^ 


=  0.  p  = 

•11 


=  2C,  = 


Ps 


00 


(26b) 


Vj  =  2C4  =  V^g. 

Assuming  Cj  and  C4,  P(p  ,  Vg)  is  calculated. 

"  *0 

(9)  The  product  P(Pgg,  Vg)  P(PJ  ,  +)  is  then  integrated  over  4*  using  Simpson's 

rule:  i.  e. , 
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*1*  *4*0 

+  4(4*1  +  I|i,  +  .  . . )  +2  {+,  +  4*4  +•  •  • )  +4'„]  ■  (27) 

(10)  P(P.  Vi")  Is  calculated  independently  from  eq  26a  and  the  result 

compared  to  eq  27.  * 

(11)  Adjustment  of  Cjand  Cels  made  until  the  values  obtained  in  eq  9  and 

10  are  equal.  Since  each  constant  catTbe  independently  adjusted  until  this  requirement 
is  met.  it  is  necessary  to  plot  C%  vs  C4  for  two  distinctly  different  sets  of  P,  and 
values.  The  intersection  of  these  two  curves  is  assumed,  without  proof  of  ulkiqueness, 
to  be  the  solution,  i.e. ,  and  Vgg. 

(12)  Since  we  are  dealing  with  saltation  only,  and  assuming  this  phenomenon 
to  be  most  responsible  for  small-scale  "drifting"  (no  large  surface  obstructions),  it 
seems  reasonable  to  assume  that  any  event  occurring  on  the  field  will  be  duplicated  a 
distance  downstream  corresponding  to  the  "most  probable"  trajectory.  This  trajectory 
is  determined  by  Pg..  and  Vgg,  and  its  "range"  may  be  interpreted  as  the  propagation 
wave  length.  It  is  calculated  with  eq  9. 

(13)  The  vertical  particle  distribution  in  the  saltation  layer  may  now  be 

obtained: 

a)  the  height  above  the  surface  is  chosen,  as  are  several  evenly 
spaced  distances,  AXg,  upstream  of  the  location  in  question. 

b)  For  various  tg's  (corresponding  to  the  position  of  a  particle  at]^g, 
Axg)  V(0)  and  ic(0)  are  calculate<ffrom  eq  5  and  9.  respectively;  -^g  and  xg  are  fixed 
here. 

c)  Having  Pg^  and  vg  for  various  _tg '  s ,  we  plot  Pg,  vs  Vo>  repeat 
this  procedure  for  other  Axg  values. 

d)  Now,  for  each  Axg,  we  may  calculate  P(Pgg,  vX)  as  a  function  of 
Pss*  and  we  essentially  have  P(Pg.,  Axg).  describing  the  number  of  particles  emitted  at 
AXg  such  that  they  must  all  pass  tnrough  .the  point  ^g  above  the  surface.  Integration 
over  Axg  gives  the  total  number  of  particles  passing  through  ^g,  Xg;  i.e. , 

w 

yxo  max  (<•  -y- 

J  “•)  o  (“) 

0  0 

which  may  be  evaluated  by  Simpson's  rule  for  two  independent  variables.  For  nine 
known  values  this  becomes: 

~  "T^*^^**  *  *  ^  **•*  ^  **•*  ****^  *  ^®*  ^•®|‘ 

h  refers  to  interval  subdivision  sixe  (since  units  are  different  for  Pg.  and  Axg,  h  =  1 
Fs  acceptable  for  purposes  of  comparison  at  different  £g's),  and  P. .  indicates  values  of 
P(Pgg,  Axg)  at  the  station  ^  as  shown  below: 
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e)  N,(  ,  ii  calculated  for  several  relative  magnitude  of 

these  values  may  be  taken  as  the  particle  distribution  at  an  arbitrary  point  far  down¬ 
stream. 


Several  details  should  also  be  included  in  the  program  just  outlined: 

a)  In  thie  quadratic  expression  defining  a,  only  one  solution,  the  nega¬ 
tive  angle,  is  physically  possible  (because  of  the  presence  of  a  surface). 

b)  Limits  are  imposed  on  by  both  the  maximum  attainable  vertical 
(terminal)  velocity  and  the  interference  of  surrounding  field  particles  (see  eq  11),  re¬ 
membering  that  pgi  depends  on  b,  and  b^{n  ^  function  of  tp. 

c)  Because  of  the  existence  of  a  terminal  velocity  (and  secondarily 
because  of  inelastic  impacts),  there  is  a  maximum  height  to  the  saltation  layer  and  a 
maximum  velocity  attainable  by  any  particle.  This  maximum  velocity  could  be  calcu¬ 
lated,  but  in  this  study,  for  purposes  of  comparison,  it  is  assumed  at  all  times  to  be 
equal  to  the  vector  sum  of  the  free  stream  and  terminal  velocities. 

d)  From  the  vector  diagram  <  90*;  it  cannot  be  "specified"  larger. 
Numerical  example 

As  yet  an  undecided  question  is  the  possibility  of  writing  drag  as  kgVf.  Since  the 
particle  is  not  at  rest  in  any  part  of  the  trajectory  analysis,  we  may  consider  the 
Reynolds  no.  range  from  10'*  to  10'*'*>  in  which  values  of  Cq,  the  drag  coefficient,  are 
obtained  from  Hoerner  (1958)  (Cq  for  the  sphere).  And,  assuming  Cq  =  lOO/Rej,  we 
obtain  a  reasonable  approximation  to  the  experimental  variation.  TRen, 


kg 


JP  c 


D^a  V,  = 


P<r. 


For  purposes  of  future  experimental  verification,  the  following  set  of  physical 
data  was  assumed: 

1.  particle  density  s  100  Ib/ft*  (commercial  glass  spheres) 

2.  particle  coefficient  of  restitution  =  0.  94 

3.  particle  diameter  =  0.2  mm. 

4.  boundary  layer:  linear  variation  to  15  ft/sec  in  6  in.  ;  i.  e.  , 
kg  s  30  (standard  atmosphere  assumed). 

A  trajectory  was  traced  out  (Fig.  9)  using  possible  values  for  Pg  and  Vg.  This 
path  appears  quite  reasonable  when  compared  to  photographs  obtained  by  Bagnold  (1941) 
for  sand. 

The  analysis  outlined  above  (p.  2-5)  gave  the  following  results: 

•)  Psoo  ”  ^00  =  8.6  ft/sec;  P(Pgg>v^)  given  in  Figure  10. 

b)  The  "ripple  wave  length"  was  then  calculated  to  be  1.6  in. 

c)  The  particle  number  density,  calculated  at  four  points,  is  repre¬ 
sented  in  Figure  11  as  a  vertical  distribution  with  arbitrary 
reference  value. 

These  results  are  also  quite  reasonable  when  compared  with  tests  using  borax  of  nearly 
0. 2  mm  average  diameter  and  similar  density.  The  major  objection  to  such  a  compari¬ 
son  is  the  lack  of  sphericity  in  these  particles.  However,  it  is  possible  to  assume  an 
equivalent  statietical  diameter  for  rotating  prismatic  shapes  (during  trajectory)  equal  to 

e  ^  value.  This  question  must  be  investigated  further,  and  more  conclusive 

results  will  be  obtained  from  the  proposed  tests  with  glass  spheres. 
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A  numerical  example  each  ae 
this  aervee  to  indicate  the  plaueibiUty 
of  the  analysia  ae  carried  out  thus  far, 
and  also  points  up  the  desirability  of 
either  simplifying  the  equations  or 
programming  them  for  automatic 
computation. 

Summary  of  analysts  of  saltation 
irftenomenon 

a)  While  certain  physical 
and  matitematical  assumptions  are 
still  open  to  question,  the  quantitative 
analysia  of  ttie  saltation  layer  is  botii 
feasible  and  physically  meaningful. 

b)  The  analysis  should  be 
extended,  approximately,  to  particles 
with  non*apherical  shape  and  varia¬ 
tions  in  diameter. 

c)  The  variation  of  all  re¬ 
sults  with  changes  in  boundary  layer 
and  field  arrangement  assumptions 
should  also  be  investigated. 

d)  The  effects  of  obstruc¬ 
tions  must  be  considered  together 
with  those  of  wavy  surfaces  in  an 
effort  to  determine  the  stable  profile 
shape  of  a  snow  surface.  Many 
interesting  confirmatory  experiments 
are  possible  in  this  phase  of  the  study. 

e)  Results  of  this  and  sub¬ 
sequent  analyses  should  ultimately  be 
applied  to  the  detailed  formulation  of 
modeling  criteria.  Thus,  for  example,  if  it  is  particle  distribution  and  wave  configura¬ 
tion  that  must  be  duplicated  in  scale  model  experiments,  it  would  be  possible  to  deter¬ 
mine  the  required  scale  factors  to  achieve  this.  Scale  experiments  then  would  afford 
the  best  means  of  time  studies  of  large-scale  drifting  and  accumulation. 

development  of  scale  factors 

Scale  factors  based  on  dimensional  analysis 

The  complex  nature  of  snow-drifting  phenomena  makes  the  determination  of 
modeling  criteria  especially  difficult.  Conflicting  requirements  are  encountered  and  it 
must  be  decided  which  may  be  neglected  without  serious  loss  of  accuracy.  Those 
criteria  which  appear  most  important  at  this  time  are  presented  in  the  following  section. 
Further  experience  with  snow-drift  modeling  will  undoubtedly  lead  to  improvements  in 
modeling  criteria. 

Scale  factors  for  establishing  modeling  criteria  will  be  developed  by  two  methods. 
Dimensional  analysis  will  be  Uie  primary  basis  for  the  first  mediod,  while  the  equations 
of  motion  of  a  particle  will  be  used  in  the  second.  The  latter  covers  part  of  the  scale 
factors  developed  by  the  first  method. 

The  development  of  the  scale  factors  will  be  carried  out  in  steps  based  on  various 
phases  of  snow-drifting  phenomena.  Most  of  the  scale  factors  will  be  developed  by  di¬ 
mensional  analysis  but  some  will  be  based  on  characteristics  of  the  physical  phenomena 
involved. 


numcLi  ocssiTv 

Figure  11.  Particle  distribution. 
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At  this  stage  of  development,  it  appears  advisable  to  concentrate  on  modeling  of 
geometric  properties.  It  is  certainly  necessary  to  preserve  geometric  configuration 
of  snow  drifts  and  for  many  applications  this  may  be  the  most  important  characteristic. 
Modeling  of  other  properties  such  as  mass  transport  and  density  of  snow  in  suspension 
may  be  achieved  automatically  if  geometric  similarity  of  detailed  characteristics  of 
the  phenomena  are  preserved,  particularly  the  paths  or  trajectories  of  individual  parti¬ 
cles.  Similarity  of  paths  will  be  the  basis  for  the  following  analysis.  It  is  obvious  that 
similarity  of  paths  will  have  little  significance  unless  there  is  also  geometric  similarity 
of  the  rigid  boundaries  such  as  buildings,  snow  fences,  roadways,  etc.  It  will  be  as¬ 
sumed  that  such  similarity  exists,  although  it  could  be  included  in  the  analysis. 

Snow  drifting  is  known  to  occur  in  three  phases:  namely,  suspension,  saltation, 
and  creep.  The  scale  factors  which  are  significant  to  the  suspension  phase  will  be  con¬ 
sidered  first.  A  snow  particle  will  be  treated  as  being  in  a  path  or  trajectory  Influenced 
only  by  fluid  and  gravitational  forces.  The  following  nine  variables  are  considered  as 
significant: 

L  =  linear  reference  dimension  of  rigid  boundary  objects 
such  as  buildings  and  snow  fences 

1  =  linear  reference  dimension  of  the  paths  of  a  snow  particle 

0-  =  diameter  of  snow  particle 

m  =  mass  of  snow  particle 

Vp  =  velocity  of  snow  particle 

V  =  ambient  air  velocity  at  the  particle 

p  =  air  density 

p  =  coefficient  of  viscosity  for  air 
g  s  acceleration  due  to  gravity. 

The  following  equation  is  formed  by  employing  the  six  dimensionless  products 
which  may  be  developed  by  dimensional  analysis  from  the  above  variables. 


1/L  may  be  interpreted  as  a  measure  of  geometric  similarity  of  particle  path  to  rigid 
boundaries  as  well  as  the  ratio  of  their  linear  dimensions.  On  the  basis  of  the  above 
equation,  similarity  of  snow  particle  paths  is  preserved  if  each  of  the  five  dimension¬ 
less  products  on  the  right  side  of  the  equation  are  maintained  at  the  same  value  in  the 
model  as  in  the  prototype.  They  are  the  scale  factors. 

The  first  scale  factor  e/L  requires  the  ratio  of  sise  of  particle  to  else  of  model 
to  be  the  same  in  model  as  in  prototype.  It  could  be  broadened  to  include  similarity  of 
shape  as  well,  but  this  would  create  almost  insurmountable  experimental  problems,  and 
it  isn't  clear  at  this  stage  whether  this  refinement  is  necessary  in  view  of  other  approxi¬ 
mations  which  will  be  made.  It  may  prove  impossible  to  maintain  r/L  in  order  to  satis¬ 
fy  other  more  important  requirements. 

The  scale  factor  m/p  «*  may  be  changed  to  the  form  pu  /p  by  replacing  m/e* 
with  pp,  mass  density  of  the  particle,  since  o*  is  proportional  to  volume.  This  ratio  is 
a  measure  of  the  buoyancy  force  acting  on  the  particle.  Since  the  density  of  snow  parti¬ 
cles  is  so  much  greater  than  air  density,  the  Imoyancy  force  is  small  pompared  with 
other  forces  in  snow-drifting  phenomena.  This  scale  factor  will  be  neglected. 

The  ratio  Vp/V  requires  the  particle  velocity  to  be  scaled  in  proportion  to  fluid 
velocity.  This  scale  factor  is  significant  for  its  relation  to  the  fluid  forces  on  a  snow 
particle.  Using  the  form  conventional  in  aerodynamics,  the  resistance  force  on  a 
particle  is  given  by 
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where 


=  Cd!  n  n 

=  dimensionless  coefficient  dependent  on  the  shape  of  the 
particle  and  on  the  Reynolds  number  Re 

V,<rp 

Re  =  — jjj - ’  Reynolds  number 

>  relative  velocity  of  particle  and  surrounding  air 

S^  ~  cross  sectional  area  of  particle. 


Relative  velocity  Yx  **  vectorial  difference  between  velocities,  V  and  Vp,  in 
scalar  form 


V_  =  k, 


1  V-k4Vp 


^  -) 


(31) 


where  kj  and  are  factors  which  depend  on  the  relative  direction  of  V  and  Vp.  When  V 
and  2L  are  inThe  same  direction,  kj  =  k4  =  1.  When  the  scale  factor  Vp/V  and  geo¬ 
metric  similarity  of  particle  paths  are  preserved,  and  are  the  same  for  model  and 
prototype.  Vj,  is,  therefore,  proportional  to  Vp  andarag  may  be  expressed  in  terms  of 


g  »-D  2  p 

where  ^  is  another  constant  which  will  have  the  same  value  in  model  and  prototype. 


(32) 


The  last  scale  factor  Vpp  ir/p  is  a  form  of  Reynolds  number.  \L  may  be  replaced 
with  Yf  in  accordance  with  the  above  discussion,  and  the  Reynolds  number  takes  the 
usualform.  Insofar  as  the  suspended  particle  is  concerned,  the  drag  coefficient  Cn 
the  only  quantity  dependent  on  the  Reynolds  number  and  must  be  selected  accordingly. 
Variation  of  Cq  with  ^  for  spherical  objects  is  well  known  (see  Schlichting,  1955, 

Fig.  1. 5)  and* will  be  used  as  the  basis  for  discussion  of  snow  particles.  Cq  varies 
inversely  with  Re  for  low  values  (Stokes's  equation  range),  while  at  relatively  high  values 
Cq  is  substantially  constant.  Re  for  snow  particles  will  generally  fall  in  between  these 
extremes,  and  no  simple  function  will  express  the  dependence  of  Cq  on  ^  for  more  than 
a  limited  range.  ~ 


The  variation  of  Cq  with  Re  will  require  the  mass  density  of  the  material  used  to 
simulate  snow  particle?  to  be  dHTerent  from  that  of  prototype  snow  if  the  scale  factor  on 
particle  sise  is  to  be  maintained.  To  maintain  equal  ratio  of  particle  acceleration 
Dg/m  (due  to  fluid  forces)  to  gravitational  acceleration  g,  the  following  product  must  be 
preserved.  On  substituting  the  value  for  Dg  from  eq  32,  it  takes  the  form 


mg  8  mg  ■ 

Since  m  is  proportional  to  pps^,  the  product  may  be  written  (omitting  the  constants) 


(33) 


(34) 
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When  the  Froude  number  is  divided  out,  there  remains 


(35) 


This  is  one  form  of  the  scale  factor  which  accounts  for  the  effect  of  the  Reynolds  number 
on  particle  fluid  forces.  An  alternate  form  expresses  it  in  terms  of  particle  free-fall 
velocity  Vf.  By  equating  the  weight  of  the  particle  to  its  drag  in  terms  of  Vf,  the 
equation  Tor  Cois  found. 


C 


D 


Tvf- 


(36) 


where  is  another  constant.  On  substituting  this  value  for  Cn  into  expression  34,  an 
alternate  form  for  the  scale  factor  is  found  (omitting  constant^,  i.  e. , 


Vp/v,. 


(37) 


In  view  of  the  scale  factor  Vp/V,  Vp/Vf  may  be  replaced  by 


Vf/V.  (38) 

This  form  has  certain  advantages  in  performing  model  experiments.  or  Cn  must  be 
determined  for  the  Reynolds  number  at  which  the  particle  moves.  The  Reynmds  number 
varies  since  relative  velocity  Vf  will  vary  and  an  average  value  must  be  taken.  In  sus¬ 
pension,  a  particle  will  have  relative  velocity  equal  to  free  fall  velocity  when  it  is  not 
accelerating  and  a  greater  or  lesser  velocity  when  under  acceleration.  Thus  the  average 
relative  velocity  tends  to  be  the  same  as  free  fall  velocity,  and  the  experimental  determi¬ 
nation  of  free  fall  velocity  for  the  particles  to  be  used  as  snow  simulators  automatically 
selects  the  correct  Reynolds  number. 

While  not  specifically  included  in  the  above  analysis,  the  wind  velocity  profile 
above  the  snow  surface  must  be  considered.  Since  the  fluid  forces  depend  on  the  local 
relative  velocity  of  air  and  snow  particles,  the  velocity  field  should  be  properly  scaled. 
Two  different  cases  must  be  distinguished.  The  first  concerns  snowdrift  formation 
around  objects  exposed  to  the  air  stream.  Air-flow  patterns  around  objects  with  sharp 
edges  such  as  buildings  and  snow  fences  remain  geometrically  similar  to  relatively 
small  scales  (low  Reynolds  number)  and  the  turbulent  velocity  components  tend  to  remain 
proportional  to  the  mean  stream  velocity.  Thus  those  phases  of  snow-drifting  phenomena 
dominated  by  this  situation  are  likely  to  be  properly  scaled  with  regard  to  the  velocity 
field. 


The  second  case  deals  with  snow  movement  over  flat  surfaces  and  presents  greater 
modeling  problems.  The  extent  to  which  velocity  profiles  over  flat  surfaces  must  be 
accurately  modeled  remains  to  be  determined.  Since  most  of  the  snow  movement  takes 
place  quite  close  to  the  surface,  it  may  be  sufficient  to  model  only  the  lowest  layers. 
Analysis  of  this  problem  is  hampered  by  lack  of  data  for  the  atmosphere  as  well  as  by 
the  controlled  airstreams  of  the  wind  tunnel  laboratory.  Such  information  as  is  available 
for  flow  over  rough  surfaces  shows  some  possibility  of  proper  modeling.  Schlichting 
(1955',  eq  20,  32a)  gives  the  following  equation  for  a  velocity  profile  over  a  rough  surface 
for  completely  "rough"  flow  (to  be  discussed  later). 

V/V^  =  2.  5  log^(y/k_)  +  8.  5  (39) 

where 

V  =  mean  velocity  at  elevation  ^ 

=  friction  velocity 

=  •‘TtoIp 
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To  =  fluid  shear  etrese  at  the  surface 
y  =  elevation  or  height  above  the  surface 

=  grain  else  of  sand  particles  used  in  the  original  experiments. 

Friction  velocity  Ve  is  dependent  on  particle  sise,  but  its  evaluation  will  depend  on  how 
a  scale  model  of  tRe  atmospheric  profile  is  to  be  produced  in  a  wind  tunnel.  There  is 
only  a  limited  amount  of  data.  There  is  also  a  fundamental  question  on  whether  velocity 
profiles  produced  over  the  almost  limitless  distances  on  the  earth's  surface  can  be 
properly  scaled  in  the  severely  limited  sise  of  the  wind  tunnel  test  section.  The  boundary 
layer  (region  in  which  the  velocity  increases  from  zero  at  the  surface  to  the  free  stream 
velocity  above)  has  little  or  no  thickness  at  the  beginning  of  the  test  section  and  increases 
in  thickness  with  distance  downstream.  There  appears  to  be  no  similar  growth  of  bounda- 
ry  layer  thickness  in  the  atmosphere.  In  addition  to  thickening  of  the  boundary  layer  on  a 
finite  length  surface,  there  is  also  a  lowering  of  surface  shear  stress  with  downwind  dis¬ 
tance.  The  surface  shear  stress  is  important  to  snow  movement  near  the  surface. 
Whether  a  boundary  layer  of  constant  thickness  and  shear  stress  is  necessary  is  not  clear 
but  rapidly  changing  characteristics  are  certainly  not  desirable  where  snow  movement 
over  extended  flat  surfaces  is  involved. 

With  a  long  wind-tunnel  test  section,  experiments  may  be  performed  at  a  down¬ 
stream  location  where  the  layer  has  reached  a  desirable  height  and  where  rate  of  growth 
and  change  in  shear  stress  are  at  sufficiently  low  values.  Artificial  thickening  of  the 
boundary  layer  may  reduce  the  length  of  test  section  needed.  If  the  test  section  is  suffi¬ 
ciently  long  (with  or  without  artificial  thickening),  the  boundary  layers  on  the  floor, 
walls,  and  ceiling  will  eventually  meet  one  another  and  the  uniform  conditions  of  a  long 
tube  will  be  reached.  This  may  yield  a  suitable  constant  thickness  and  constant  shear 
stress  boundary  layer  but  this  subject  requires  further  exploration.  That  the  shear 
stress  is  constant  when  the  uniform  conditions  of  a  long  tube  are  reached  is  shown  by 
the  equation  for  shear  stress  given  in  Schlichting  (1955,  eq  20,  35). 

K  =  - i - -  (40) 

[2  log,*(R/k^)  +  1.  74)* 

where  K  is  a  coefficient  of  resistance  defined  by  the  equation 


To  =  X  p  u**/8. 


(41) 


u  is  the  mean  velocity  of  flow  through  the  tube,  R  is  the  radius  of  the  tube  and  k  is  the 
grain  size  as  introduced  above.  Since  the  maximum  velocity  will  occur  at  the  <^nter  of 
the  tube,  R  is  the  boundary  layer  thickness.  For  non-circular  shapes.  R  will  be  pro¬ 
portional  to  boundary  layer  thickness.  Preservation  of  geometric  similiTrity  of  snow 
particles  in  relation  to  boundary  layer  thickness  or  other  measures  of  velocity  profile 
dimensions  will  be  equivalent  to  keeping  R/k^  constant  between  model  and  prototype. 
This  will  yield  the  same  resistance  coefficient  as  eq  40.  In  terms  of  friction  velocity 


V*  =  V  To/p 

v^,  = 


(42) 


This  equation  has  special  significance  when  viewed  in  relation  to  eq  39  for  velocity  pro¬ 
file.  On  substituting  the  value  for  Vs  from  eq  42 


V,'a  "  -JfJf  [  2.  5  log  (y/k  J  +  8.  5] .  (43) 

Since  at  a  geometrically  similar  elevation  )r,  y/kg  will  be  the  same  for  model  and  pro¬ 
totype,  the  velocity  profile  will  be  properly  scaled. 

The  saltation  phase  of  snow  movement  involves  complex  interactions  between  snow 
particles  which  depend  on  their  elastic  or  plastic  properties.  The  rebound  of  a  particle 
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after  impact  with  one  or  more  particlea  or  a  boundary  eurface  is  an  important  charac¬ 
teristic.  The  ratio  of  velocity  of  rebound  to  velocity  of  impact  is  defined  as  the  coeffi¬ 
cient  of  restitution  e.  Since  it  is  desired  to  preserve  the  velocity  field,  the  coefficient 
of  restitution  should  be  the  same  for  model  and  prototype;  e  becomes  another  scale 
factor  in  the  modeling  of  snow  drifting.  ~ 

The  creep  phase  of  snow  drifting  involves  gravitational  and  fluid  forces.  When  a 
snow  particle  at  rest  is  lodged  among  other  snow  particles,  the  weight  of  the  particle 
produces  a  moment  which  opposes  the  overturning  moment  produced  by  air  movement. 
When  the  aerodynamic  moment  equals  the  weight  moment,  the  particle  is  on  the  thres¬ 
hold  of  movement.  The  velocity  at  which  this  occurs  is  the  threshold  velocity.  Bagnold 
(1941)  has  developed  the  following  equation  for  threshold  friction  velocity  for  sand  move¬ 
ment.  It  will  have  the  same  form  for  snow: 


^eI!. 


g  <r 


(44) 


where  is  a  constant  found  to  have  a  value  of  0.  1  for  sand  in  air  and  0.  2  for  sand  in 
water.  Its  value  is  not  known  for  snow.  If  the  ratio  pn/p  has  the  same  value  for  model 
and  prototype,  eq  44  will  have  the  same  form  as  the  Froude  number  and  will  be  pro¬ 
portional  to  air  velocity  V.  Owing  to  viscous  effects  on  particle  drag,  p^/p  will  change 
in  accordance  with  eq  35. 


Equation  44  holds  only  for  "rough"  flow,  i.  e.  ,  t’  condition  for  which  the  turbu¬ 
lent  boundary  layer  is  close  enough  to  the  surface  to  c.xvelop  the  particles.  As  the  rough¬ 
ness  Reynolds  number. 


Rer  = 


V^k.p  V^k. 


(45) 


becomes  smaller,  a  condition  is  reached  where  the  laminar  sublayer  has  a  thickness 
comparable  to  grain  size  and  eq  44  is  invalid.  A  minimum  V*t  i*  reached  and  for 
smaller  particles  begins  to  increase.  Dependence  of  pp/p  in  eq  44  and  the 

existence  of  a  minimum  Vat  laminar  sublayer  may  make  the  model  Vat  higher  than 

that  which  satisfies  the  froude  number.  Thus  atmospheric  wind  velocities^  close  to 
atmospheric  threshold  velocity  may  be  difficult  to  model.  This  situation  is  improved  by 
the  fact  that  once  snow  movement  has  started,  it  will  continue  at  a  velocity  lower  than, 
threshold  velocity,  called  the  impact  threshold  velocity.  Feeding  the  simulating  snow 
from  the  ceiling  of  the  test  section  in  a  wind  tunnel  experiment  will  reduce  the  threshold 
velocity.  The  presence  of  objects  causes  local  regions  of  higher  velocity  where  snow 
movement  may  start  at  a  lower  threshold. 


Atmospheric  threshold  velocity  for  snow  movement  is  sometimes  given  as  10  mph. 
Assuming  this  velocity  occurs  at  an  elevation  of  6  ft  and  the  snow  particle  diameter  is 
0.  1  cm,  the  atmospheric  threshold  shear  velocity  may  be  calculated  with  eq  39.  Vet 
found  to  be  16.  4  cm/sec.  The  roughness  Reynolds  number  corresponding  to  these 
conditions  is  approximately  12,  assuming  grain  size  is  the  same  as  particle  diameter. 
Bagnold  (1941)  reports  that  the  flow  around  sand  partl»es  ceases  to  be  "rough"  when 
Re^  is  less  than  3.  5.  Schlichting  (1955)  shows  the  flow  to  be  "smooth"  when  is  less 
than  5  and  "rough"  when  greater  than  70,  with  a  transitional  region  in  between.  His 
data  are  based  on  experiments  in  which  particles  are  not  free  to  move.  While  there  may 
be  differences  in  the  results  due  to  differences  in  method  of  measuring  particle  size, 
the  large  difference  between  Schlichting's  experiments  and  the  others  is  that  the  surface 
of  loose  snow  or  sand  may  have  been  pitted  by  prior  saltation  leaving  craters  larger  than 
particle  diameter.  This  would  make  the  "effective"  particle  size  larger.  Regardless  of 
these  differences,  Re^  for  snow  in  the  prototype  is  not  much  greater  than  that  for  mini¬ 
mum  Vet,  il^us  lea^^g  little  room  for  reducing  model  scale  without  limiting  velocity 
rangeoecause  of  increasing  Vet*  U  the  lower  part  of  the  velocity  range  becomes  impor¬ 
tant,  it  may  be  advisable  to  consider  abandoning  the  requirement  of  geometric  similarity 
for  particle  size.  This  will  obviously  require  change  in  particle  density  and  may  give 
greater  freedom  in  the  selection  of  material  for  snow  simulation. 
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In  aummary,  the  following  scale  factors  are  considered  to  be  significant  at  this 
stage  of  development: 


V* 

0-  _p 

I.- 


’'f 

■7- 


e. 


In  application,  the  second  and  third  factors  may  be  combined  to  form  a  Froude  number 
based  on  fluid  velocity  V*/g(r.  The  density  of  the  model  particle  must  be  such  as  to 
satisfy  both  (t/L  and  Vf/V.  Consideration  should  be  given  to  proper  scaling  of  the  wind 
velocity  profile,  especially  near  the  surface.  The  discussions  of  velocity  profiles  in 
this  paper  are  baaed  on  characteristics  obtained  when  there  is  no  movement  of  the  parti¬ 
cles.  When  the  snow  particles  are  in  motion,  the  velocity  profile  is  different.  Whether 
the  above  modeling  criteria  are  adequate  for  this  situation  is  a  subject  for  further  explo¬ 
ration. 


Scale  factors  baaed  on  theoretical  analysis  of  saltation  phenomena 

The  theoretical  analysis  provides  a  basis  for  determining  scale  factors  which 
apply  to  saltation  phenomena.  The  assumptions  of  that  analysis  will  obviously  apply 
here.  The  surface  configuration  is  that  of  a  flat  uniform  field  of  particles.  The  follow¬ 
ing  analysis  begins  with  the  differential  equations  of  motion  which  apply  to  both  model 
and  prototype.  The  condition  for  similarity  of  trajectories  is  applied  and  the  scale 
factors  needed  to  obtain  this  similarity  are  found. 


J  ,  ^  t  ,  .  0  (46) 

(«) 

If  (  )_  denotes  model  variables  and  we  introduce  the  scale  factors: 

'nfi 

length:  L,  =  (48) 

velocity:  U.  =  -r-  (49) 

•  L 

particle  density:  D^  = 


we  may  calculate  the  factor  for  acceleration  quite  simply,  assuming  a  steady  velocity 
field: 


•  m 


or. 
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The  equations  of  motion  are  now  written  for  the  model,  and  eq  48,  49,  50  and  51 

substituted,  giving: 


(since  ko  ~  tr) 


k,L.^s 

mD^L* 


■  mD.L* 


iL,y 


where  K  =  (52) 

^1 


However,  had  we  selected  a  different  scale  factor  than  Lg  for  particle  siae  (leaving  Lg 
for  dimensions  that  are  large,  compared  to  particle  diameter),  we  would  have  obtained. 


(53) 

(54) 


The  procedure  is  certainly  allowable  when  (a)  particle  diameter  is  much  smaller  than 
the  critical  dimensions  of  large  objects  placed  on,  or  restricting,  the  field,  and  (b) 
particles  do  not  influence  each  other:  i.  e. ,  collisions  between  airborne  particles  are 
unimportant.  Then,  for  a  model  particle  to  behave  like  a  full-scale  saltating  particle, 
eq  53  and  54  must  reduce  to  eq  46  and  47  respectively.  It  is  obvious,  then,  that 


and 


KLs 

15:^ 


1 


are  required  of  the  model. 

Before  proceeding,  we  note  that  the  factor  k  is  certainly  not  an  independent 
parameter.  Since 


m 


>^i 


then, 


(55) 

(56) 

(57) 


(58) 
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Therefore,  condition  (12)  reduces  to: 


U  L 

•  s 


U. 


L  D  r*  =  P  o>^  to  eq  56. 


And,  in  satisfying  eq  55  and  56,  we  get 


(59) 


[Written  as  V|n/lm8  =  V*/lg,  eq  59  is  recognized  as  the  statement  that  the  Froude 
number  is  constant.  J 


and 


From  eq  58,  we  get: 


(60) 


(61) 


Our  modeling  must  now  be  subjected  to  a  test  for  self-consistency.  That  is,  we  have 
required  that: 


-U 

IT,- 


and 


(mg), 


m 


mg 


But,  our  basic  eq  46  and  47  have  been  derived  using  the  assumption  that 
Drag  =  k,u^, 
or 

kjVj  =  mg 


J^fm  .  <»"g>m  Xk.  ,  _  f  ^  , 
Vf  mg  \\<om  /  ■  IrngT  \<r^/ 

which  implies  that 

U,  =  D,  E*/E  =  D,  r* 

i.  e.  , 

^  “• 

But  this  condition  has  already  been  met  by  eq  56. 
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Therefore,  the  scale  factors  are  consistent  with  the  physical  assumptions. 

The  model  is  selected  as  follows: 

(1)  A  convenient  ratio  L _ /L  is  chosen. 

(2)  Free  stream  and  particle  velocities  must  now  conform  to  the^atio 


the  boundary  layer  being  adjusted  so  that 


(3)  The  particle  density  and  size  must  be  selected  so  that 


There  is  some  freedom  in  selecting  these  properties,  but  an  upper  limit  is  imposed  on 
density  by  the  requirement  that  the  particle  possess  a  reasonable  threshold  velocity 
(conforming  to  the  magnitude  of  other  velocities  in  the  problem).  The  lower  limit  on 
density  and  upper  limit  on  size  are  obvious. 


Our  modeling  scheme  is  now  virtually  complete.  We  need  only  consider  the 
boundary  conditions;  that  is,  we  must  require  that  termination  of  one  trajectory  at 


angle  i)<,  and  velocity  yi  leads  to  the  initiation  of  another  trajectory  at  angle  Po  veloci¬ 
ty  'I'lm  =  't'l*  Pom  =  Po 


(from  the  condition  tan 


and  Vjin  =  ^i>  ^om  =  Uvg'.  The  exact  requirement  on  the  ratio  of  coefficients  of 

restitution,  em/e,  ao  that  these  conditions  are  met  is  determined  by  inspection  of  the 
interaction  equations.  These  equations  show  that  the  variables  in  question  are  inde¬ 
pendent  of  particle  velocity  and  size;  hence  the  model  interaction  will  satisfy  the  above 
conditions  only  when  em/a  =  f- 


Selection  of  model  scale  and  snow  simulator  properties 


The  criteria  developed  in  the  preceding  sections  are  used  to  select  the  require¬ 
ments  to  be  applied  in  designing  the  scale  model  experiments.  At  the  earliest  stage  of 
the  project,  before  these  criteria  were  developed,  the  search  for  a  suitable  snow  simu¬ 
lator  was  based  on  finding  a  material  which  had  properties  similar  to  real  snow.  Light¬ 
weight  materials  seemed  to  satisfy  this  requirement  and  the  early  tests  in  the  10  x  20  in. 
wind  tunnel  used  materials  of  this  type. 

In  order  to  select  the  most  effective  linear  scale  for  the  physical  objects  to  be 
used  in  the  3^  by  7  ft  wind  tunnel,  consideration  had  to  be  given  to  the  size  of  the  object 
which  can  be  placed  in  the  tunnel  without  being  seriously  affected  by  the  wall  and  ceiling 
boundary  conditions.  Based  partly  on  the  fact  that  the  snow  drifting  studies  are  cdn- 
cerned  with  surface  flow,  it  was  felt  that  objects  as  much  as  2  ft  high  could  be  used 
without  any  serious  errors  being  introduced.  Consideration  had  to  be  given  also  to  how 
small  the  models  could  be  and  still  achieve  meaningful  results.  If  the  models  and  model 
snow  are  too  small,  then  the  boundary  effects  and  other  random  effects  might  be  of  the 
same  order  of  magnitude  as  the  desired  drifting  effects.  Another  factor  considered  was 
the  available  wind  velocity  in  the  wind  tunnel.  This  limits  the  size  of  the  model  snow 
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particles  aince)  aa  shown  later,  the  fluid  threshold  velocity  for  very  small  particles 
is  an  inverse  function  of  the  particle  siae.  After  review  of  the  problem,  it  was  decided 
to  use  an  initial  linear  scale  of  1/10.  These  conditions  would  allow  the  modeling  of 
objects  which  in  reality  had  heights  of  as  much  aa  20  ft  and  widths  of  up  to  approximately 
40  to  50  ft.  Since  natural  blowing  snow  on  the  Greenland  Ice  Cap  has  been  reported 
(Gerdel,  personal  communication)  to  have  a  mean  grain  diameter  of  1  mm,  the  simu- 
lating-snow particles  must  have  a  mean  diameter  of  0. 1  mm. 

The  air  velocity  requirement  is  determined  from  the  Froude  number  aa  follows: 


gL  -  gL^ 
or  for  constant  g. 


Thus  for  a  1/10  linear  scale  the  velocity  must  be  modeled  •/I/To. 

The  free  settling  velocity  of  the  snow- simulator  particles  must  also  satisfy  the 
Froude  number  and  must,  therefore,  have  a  value  ■s/l/lO  of  that  of  real  snow  particles. 
The  density  of  the  simulator  should  satisfy  the  particle  diameter  and  free  settling  veloc¬ 
ity  requirement  simultaneously.  The  coefficient  of  restitution  of  the  simulator  should 
be  the  same  as  that  of  real  snow. 

If  the  velocity  profile  were  modeled  as  discussed  earlier,  Ve/V,  the  ratio  of 
friction  velocity  (a  measure  of  surface  shear  stress)  to  airstream  velocity  (at  geo¬ 
metrically  similar  elevations),  would  remain  constant  as  shown  by  eq  39.  For  proto¬ 
type  elevation  y  of  6  ft  and  snow  particle  diameter  of  0. 1  cm  as  used  earlier,  eq  39 
gives  Ve/V  -  .^37.  It  is  of  interest  to  compare  this  with  a  value  estimated  for  the 
model  experiments.  Velocity  profile  measurements  were  not  made  in  these  experi¬ 
ments  but  an  approximated  value  may  be  obtained  with  the  equation  for  surface  friction 
coefficient  over  a  rough  plate  (Schlichting,  1955,  eq.  21.41): 

C/  »  1 2.  87  +  1 . 58  In  (x/k^)J  »  (62) 

where  ^  is  the  local  friction  coefficient  and  x  the  distance  downstream  of  the  leading 
edge  ofthe  plate.  C{  is  related  to  Ve,  by  the  equation, 

V*  = 

where  y«  is  the  free  stream  velocity  above  the  boundary  layer.  The  wind  tunnel  experi¬ 
ment  does  not  meet  the  conditions  of  eq  62  because  there  is  a  section  of  smooth  floor 
upstream  of  the  snow  surface,  but  an  approximate  value  may  be  obtained.  C{  changes 
slowly  with  x.  Assuming  an  effective  value  of  15  ft  for  x  and  0.  01  cm  for  particle  siae, 
eq  62  and  63  yield  Ve/V  =  0.  039.  This  is  comparable  t^ the  value  calculated  for  the 
atmosphere.  For  a  particle  diameter  of  0.  02  cm,  Ve/V  becomes  0.  042.  There  arises 
the  question  of  whether  geometric  similarity  of  profiles  will  be  preserved.  With  equal 
values  of  Ve/V,  geometric  similarity  is  implied  in  eq  39  but  there  is  the  limitation  of 
boundary-layer  depth.  No  equation  for  boundary-layer  depth  appears  to  be  available  for 
the  rough  plate  condition  of  eq  62.  The  depth  is  estimated  to  be  on  the  order  of  |  ft  for 
X  =  15  ft.  The  model  experiments  may,  therefore,  be  expected  to  give  approximate 
modeling  of  the  velocity  profile  for  a  depth  corresponding  to  5  ft  in  the  atmosphere  at  x  & 
15  ft.  Above  this  level  the  model  air  flow  is  at  constant  velocity,  while  in  the  atmos¬ 
phere  the  velocity  continues  to  increase  to  a  much  greater  elevation. 

Modeling  of  fluid  threshold  velocity  for  snow  movement  presents  the  problem  dis¬ 
cussed  earlier  with  regard  to  lower  limit  on  rough  flow.  A  minimum  roughness  Reynolds 
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number  of  5  will  be  taken  as  given  by  Schlichtlng  (1955)  for  the  beginning  of  smooth 
flow.  For  particle  diameter  of  0.  01  cm.  Vet  ^  cm/ sec  (eq  45),  but  for  0.  02  cm; 

Vet  ~  cm/sec.  The  free-stream  threshold  velocity^  fpr  Ve/V  *  0.  042  and  0.  02  cm 
?iam  is  19  mph.  The  equivalent  prototype  speed  woulifoe  very  high.  If  this  were 
truly  the  threshold  speed,  flow  over  flat  snow  surfaces  could  only  be  modeled  for  rather 
high  prototype  wind  speeds.  The  wind-tunnel  experiments  showed  that  snow  movement 
began  at  free  stream  speeds  on  the  order  of  1 1  mph  in  the  flat  region  upstream  of 
misled  objects  when  no  snow  was  falling  (fluid  threshold  velocity)  and  at  a  lower  speed 
with  falling  snow  (impact  threshold  velocity).  Irregularities  in  the  snow  surface  proba¬ 
bly  played  an  important  part  by  raising  the  effective  grain  sise. 

INVESTIGATION  AND  SELECTION  OF  MATERIAL  TO  BE 

USED  AS  SNOW  SIMULATORS  IN  MODEL  EXPERIMENTS 

Techniques  and  equipment 

Measurement  of  free-fall  velocity. 

Large  drop  chamber  with  still  camera  and  strobadume!  The  free-fall  velocity 
is  one  of  the  niore  important  characteristics  of  particles.  Equipment  and  techniques 
were  developed  to  measure  this  property.  In  the  first  method  tried,  a  test  chamber 
with  a  cross-sectional  area  of  10  by  20  in.  and  a  height  of  5  ft  was  used.  This  chamber 
had  Plexiglas  on  two  adjacent  sides  and\i(as  equipped  with  a  vertical  scale.  A  still 
camera  and  a  flashing  light  (strobalume)  ^re  placed  near  the  bottom  of  the  chamber  so 
that  the  light  illuminated  falling  particles  at  an  angle  with  respect  to  the  camera.  The 
flashing  rate  of  the  light  was  controlled  by  a  high-speed  switch.  With  the  shutter  of  the 
camera  open  and  the  light  flashing,  particles  were  dropped  in  the  test  chamber.  This 
produced  a  series  of  successive  images  as  the  particles  fell  past  the  camera.  Since  the 
hashing  rate  of  the  light  was  known  and  the  distance  traveled  was  measured  on  the  verti¬ 
cal  scale,  the  free-fall  velocity  td  the  particles  could  be  calculated.  In  order  to  insure 
that  free-fall  velocity  had  been  reached,  the  fall  path  of  the  particle  was  about  5  ft  long 
before  it  entered  the  camera's  view.  The  major  deficiency  of  this  method  is  that  one 
cannot  be  certain  that  two  successive  images  on  the  photographic  plate  were  made  by  the 
same  falling  particle. 

Two-camera  technique;  Two  still  cameras  and  a  time-delay  mechanism  were 
used  with  tbe  drop  chamber  just  described.  The  two  cameras  were  placed  at  different 
heights  outside  the  chamber  with  the  lower  camera  triggered  a  split  second  after  the 
upper  one.  A  millisecond  timer  measured  the  delay  between  the  two  exposures.  The 
delay  was  varied  until  it  was  possible  to  observe  the  same  particle  in  the  photos  taken 
by  the  two  cameras.  Although  the  basic  procedure  was  sound,  this  method  was  not  used 
extensively  after  a  high-speed  motion  picture  camera  became  available. 

Drop  chamber  and  high-speed  motion  picture  camera;  A  drop  chamber  was 
constructed  which  provided  back  ligating  of  tbe  /ailing  particles  by  a  1 50  w  light  source 
placed  behind  a  frosted  glass  plate  (Fig.  12).  The  falling  particles  were  photographed 
by  a  high-speed  motion  picture  camera,  *  at  approximately  600  frames  per  sec.  As 
used  in  previous  techniques,  a  strobalume  illuminated  the  camera's  field  at  a  rate  of 
50  flashes  per  sec.  In  this  way  the  flashing  light  would  appear  on  the  film  about  once 
every  12  frames.  This  method  yielded  very  satisfactory  results  and  was  used  to  deter¬ 
mine  the  fall  velocity  of  snow  grains.  A  later  refinement  to  this  method  was  made  by 
using  a  neon  light  bulb  in  the  field  of  view.  The  film  speed  was  fast  enough  to  detect  the 
on  and  off  periods  of  the  light  bulb  vdien  an  ordinary  60-cyclo  a-c  current  was  used  for 
illumination.  This  provided  a  much  better  timing  mechanism  than  had  been  used  previ¬ 
ously. 

Density  measurement.  Since  the  density  of  a  small  particle  is  an  important  fac¬ 
tor  in  its  terminal  fall  velocity,  attempts  were  made  to  develop  a  method  for  measuring 
the  densities  of  the  substances  which  might  be  suitable  for  a  model  snow.  It  was  found 
that  for  many  non-crystalline  or  non-uniform  substances  like  cork,  wood,  fly  ash,  etc. 
the  densities  reported  in  the  literature  varied  greatly  depending  on  the  conditions  under 
which  they  were  measured.  A  simple  experiment  showed  that  cork  particles  which  had 


*Bell  and  Howell  Traid  70. 


Figure  13.  Terminal  velocitiee  of  spheri¬ 
cal  particles  of  different  density  settting  in 
Figure  12.  Drop  chamber  showing  the  Bell  air  at  70F  under  the  action  of  gravity, 

and  Howell  high-speed  camera  in  position.  (Extracted  from  Perry,  1990,  p.  1021). 


been  reduced  to  powder  would  not  float  on  water  (neglecting  surface  tension  effects) 
although  the  density  of  cork  is  reported  to  be  0.  207  by  Hodgman  (1952).  Three  encroach¬ 
es  to  the  density  measurement  problem  were  undertaken. 

Weight  and  volume  technique;  The  first  approach  was  an  attempt  to  refine  the 
standard  weight-per-unU-volume  technique.  Ten -milliliter  samples  of  several  sixes  of 
test  materials  were  vibrated  until  they  appeared  to  be  well  settled  (i.  e. ,  little  measur¬ 
able  reduction  of  volume  seemed  obtainable).  Using  the  approach  described  by  Dalle 
Valle  (1948)  to  estimate  the  percentage  of  voids,  the  volume  of  the  test  solids  was  calcu¬ 
lated.  The  samples  were  then  weighed  on  a  microbalance  and  the  density  was  computed 
by  the  elementary  relationship: 

density  =  mass/volume. 

The  uncertainties  of  this  method  are  quite  apparent.  The  equations  mentioned  in  Dalle 
Valle's  method  involve  certain  assumptions  of  sphericity  and  "closeness"  of  packing. 

Since  most  of  the  test  materials  contain  unspherical  particles,  some  errors  in  die  results 
are  to  be  expected.  When  materials  of  known  densities  were  used  in  this  method  the  re¬ 
sults  appeared  to  be  too  low,  probably  due  to  incomplete  packing  of  the  sample. 

Air  suspension  technique:  Anodter  method  used  to  calculate  densities  involved 
mounting  a  large  ^ss  tube  in  a^rtical  position.  Test  particles  dropped  in  the  tube 
were  held  in  mid-ti^  by  forcing  a  stream  of  air  throu|^  die  tube.  The  velocity  of  the 
air  was  measured  by  a  ffowrator.  The  velocity  required  to  exactly  balance  the  particles 
represents  a  measure  of  the  terminal  velocity  of  the  particles.  The  densities  were  then 
calculated  using  Stokes's  law. 

Vft  s  terminal  settling  velocity  of  particle 

pp  s  density  of  particle 
p  =  density  of  fluid 

r  =  radius  of  particle  (assumed  to  be  spherical) 

g  3  gravity 

|i  s  viscosity  of  fluid. 
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An  error  ie  introduced  in  the  epplieetion  of  this  equation  to  the  problem  since  the  actu¬ 
al  curves  become  non-linear  for  particle  diameters  greater  than  50  to  75p  (Fig.  13). 

Thus  in  the  intermediate  range  of  diameters  between  SOp  and  a  few  thousand  microns 
neither  Stokes's  law  nor  the  ordinary  laws  of  free-falling  bodies  are  completely  valid. 

For  this  reason  it  was  decided  to  measure  the  terminal  fall  velocity  of  the  materials 
selected  for  possible  use  as  snow  simulators. 

Two-liquid  technique;  From  among  the  samples,  pairs  of  particles  were 
selected  so  that  each  pair  was  closely  matched  in  sise  and  shape.  One  of  each  pair 
was  dropped  into  a  column  of  distilled  water  and  its  time  of  fall  (or  rise  for  light  ma¬ 
terials)  was  measured  through  a  known  distance.  The  other  particle  was  dropped  into 
a  column  of  another  liquid  of  known  density  (carbon  tetrachloride  was  used  in  these 
experiments),  and  its  time  of  fall  was  measured  through  an  equal  distance.  The  ratio 
of  the  two  time  periods  was  the  value  sought.  A  graph  was  devel<^d  by  measuring 
this  ratio  for  a  number  of  substances  of  known  density.  This  meAod  was  abandoned 
because  of  a  revised  theoretical  approach  to  the  entire  snow  material  selection  problem, 
which  is  discussed  in  another  section. 

Free-fall  velocity  of  snow  particles 

A  literature  search  for  the  mean  terminal  fall  velocity  of  natural  snow  did  not 
yield  any  conclusive  value  for  the  type  of  snow  usually  found  on  the  Greenland  Ice  Cap. 

It  was  therefore  decided  to  set  up  an  experiment  to  determine  this  value  for  natural 
blowing  snow  as  well  as  for  all  the  proposed  model  materials.  It  was  felt  that  a  photo¬ 
graphic  technique  would  yield  the  best  results  in  this  type  of  study.  Several  arrange¬ 
ments  involving  a  drop  chamber  and  photographic  equipment  were  tried  before  satis¬ 
factory  results  were  obtained  using  a  high-speed  motion  picture  camera. 

The  experimental  phase  of  this  work  was  performed  by  New  York  University 
scientists  at  SIPRE's  cold  laboratory  in  Wilmette,  Illinois.  A  sample  of  snow  from  the 
Greenland  Ice  Cap  was  used.  The  test  was  carried  out  in  the  22F  cold  room  using  a 
Fastax  camera.  The  basic  procedure  is  described  above  although  some  modifications 
in  camera  speed  and  setting  were  made.  Several  rolls  of  high-speed  motion  pictures 
were  taken  of  falling  snow  particles  between  1. 0  and  1.4  mm  in  diameter.  On  the  re¬ 
sulting  films  20  representative  particles  were  tracked  for  1/40  see  (50  frames).  There 
appeared  to  be  great  uniformity  in  the  sise  and  velocity  of  the  snow  grains.  The  results 
are  shown  in  Table  I.  The  mean  vertical  velocity  computed  from  these  measurements 
was  209  cm/ se'c . 

Physical  properties  of  particulates  to  be  used  as  snow  simulators 

Low-density  materials.  In  the  early  stage  of  the  investigations  a  number  of  low- 
density  materials  were  selected  for  testing  as  potential  snow  simulating  materials  and 
also  to  determine,  if  possible,  whether  the  test  procedures  were  adequate.  The  results 
of  these  tests  on  a  group  of  materials  are  presented  in  Table  II.  None  of  the  materials 
tested  were  considered  to  be  suitable. 

During  the  course  of  these  tests,  it  became  apparent  that,  for  some  of  the  non- 
crystalline  materials,  there  was  a  relationship  between  particle  sise  and  density  and 
that  this  phenomenon  should  be  investigated  before  further  consideration  be  given  to  such 
materials.  Physical  inspection  showed  undesirable  characteristics  such  as  the  inabili¬ 
ty  to  remain  free  flowing  and  extreme  variability  in  physical  properties.  The  remainder 
(k  the  materials  were  tested  for  free-falling  velocity  and  density  by  one  or  more  of  the 
methods  described  previously  (Table  II).  It  was  concluded  that  further  investigation 
should  be  made  on  these  and  other  possible  materials. 

A  RoTap  Siever  and  a  set  of  U.  S.  Standard  Sieves  were  used  to  prepare  test 
samples  of  cork  particles.  A  sample  of  ground  cork  was  sieved  until  there  was  no  ap¬ 
preciable  change  in  the  quantity  of  particles  remaining  on  each  of  the  sieves.  The  sam¬ 
ples  were  then  removed  from  the  sieves  and  each  sample  was  then  resieved  briefly 
through  the  next  larger  sieve  in  the  series.  This  double  sieving  was  done  on  the  assump¬ 
tion  that  the  non-spherical  particles  may  go  through  the  sieve  endwise  over  a  long  siev¬ 
ing  period  while  a  short  sieving  period  will  retain  some  of  the  larger  ones.  The  particles 
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Table  I.  Fall  of  20  enow  grains  in  1/40  tec. 


Particle  no. 

Vertical  distance 
(cm) 

Particle  no. 

Vertical  distance 
(cm) 

1 

4.36 

11 

6.43 

2 

5.23 

12 

5.67 

3 

4.17 

13 

4.56 

4 

5.76 

14 

5.56 

5 

5.80 

15 

5.88 

6 

5.44 

16 

5.23 

7 

4.  72 

17 

5.28 

8 

5.04 

18 

6.20 

9 

5.  72 

19 

5.  56 

10 

4.36 

20 

4.88 

Table  II.  Physical  properties  of  materials 
originally  proposed  as  snow  simulators. 


Bulk  density 

Free  |all 

Exp  density 

Sise 

velocity 

Name 

(  lb/ft»)* 

(g/cm>)  t 

(mm)** 

(cm/sec)tt 

Vermieulite 

10 

0.16 

0.2 

18.  5 

Styrofoam 

1.8 

0.03 

4.0 

28.7 

Balsawood  sawdust 

7.49-12.49 

0.12 

0.3 

26.6 

Mica 

165 

2.7 

0.1 

Polystyrene 

2 

0.03 

Flocking  (cellusuede 
rayon)  ‘ 

12 

0.20  1. 

0  (length) 

13.3 

Cork 

12.92 

0. 38 

0.4 

15.  5 

Foam  rubber 

6 

0.1 

Fly  ash 

95  (variable) 

0.5 

0.1 

*  Taken  from  literature, 
t  Measured  by  weight  and  volume  technique. 

**  Moan  siae  as  determined  by  a  set  of  standiard  sieves, 
tt  Measured  by  the  still  camera  and  strobalume  method. 
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which  remained  on  the  larger  sieve  after  the  second  sieving  were  assumed  to  have  ap* 
proximately  die  dimension  of  the  larger  sieve  and  these  were  used  in  the  density 
testing.  The  samples  were  tested  by  both  the  weight-and-volume  method  and  the  air* 
suspension  method.  The  results  are  shown  in  Table  ni  and  Figure  14.  The  approach 
of  the  density  of  the  smaller  particles  toward  higher  values  agrees  with  the  observation 
of  cork  powdsr  in  water  mentioned  previously. 

High»denslty  materials.  The  following  characteristics  were  taken  as  representa¬ 
tive  of  natural  blowing-snow  particles  on  the  Greenland  Ice  Cap. 

Diameter  1  cm  (Gerdel,  personal  communication) 

Free-fall  velocity  =  209  cm/sec  (from  experiments  reported  above) 
Coefficient  of  restitution  =  0. 33  (from  experiments  reported  below) 

Applying  the  modeling  criteria  to  a  linear  scale  of  1/10,  the  following  desired 
properties  cd  the  snow  simulator  are  found. 

Diameter  =  0.  01  cm 
Free-fall  velocity  =  66  cm/sec 
Coefficient  of  restitution  =  0.33. 

In  order  to  satisfy  the  requirement  on  size  and  free-fall  velocity,  the  simulator 
material  must  have  a  density  greater  dian  that  of  snow  particles.  A  guide  to  the  approxi¬ 
mate  density  required,  which  is  needed  in  the  preliminary  selection  ^  materials,  can 
be  obtained  by  applying  the  Stokes's  law  form  o£  free-fall  velocity  equation  for  spheres. 
This  equation  gives  a  density  of  2.  1  g/cm*  for  the  above  simulator  properties.  Ma¬ 
terials  having  density  near  ^is  value  were  selected. 

Free-fall  velocity;  The  use  of  the  high-speed  motion  picture  technique  des¬ 
cribed  previously  made  it  possible  to  obtain  reliable  values  for  free-fall  velocity.  A 
series  of  tests  showed  that  the  camera  speed  was  sufficiently  constant  so  that  no  timing 
device  was  necessary.  The  limits  of  variability  of  the  camera  were  found  to  be  from 
211  to  221  frames  per  second  after  the  initial  acceleration  period  of  4/5  sec.  Since 
these  deviations  produced  a  maximum  error  of  only  2.  3%,  the  assumption  that  the 
camera  ran  at  an  average  of  216  frames  per  second  was  considered  satisfactory. 


Figure  14.  Relationship  between  particle 
diameter  and  particle  density  for  cork. 


Figure  15.  Relationship  between  particle 
diameter  and  free-faU  velocity  for  larger 
sisad  particles. 
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Table  111.  Density  tests  on  cork  particles. 


U.  S.  Standard 

Weight  and  vol. 
Siae  method 

Air  suspension 
metttod 

Mean 

Sieve  no. 

(g/cm*) 

{g/cm*i 

(g/cm*) 

170 

powder 

88 

0.87 

0. 93 

1.0(«^ 

0.90 

140 

105 

0.83 

0.  93 

0.  88 

50 

297 

0.51 

0.  59 

0. 55 

35 

500 

0.35 

0.39 

0.37 

20 

840 

0.30 

0.38 

0.34 

18 

1000 

0.20 

0.  32 

0.28 

bulk 

*  Powdered  cork  was  observed 
t  Reported  by  Hodgman  (1952). 

to  sink  in  water. 

0.207  (t) 

Table  IV.  Terminal  velocity  tests  using  high-speed  motion  picture  camera. 


Chemical 

Siae 

Density* 

Terminal  velocity  i 

Material 

formula 

(mm) 

(g/ cm») 

(cm/sec)  1 

Alum 

KAl>(SOlt),‘  12H,0 

0.6 

1.75 

350  j 

0.  1 

33  ! 

Borax 

Na,  BeOf  lOHaO 

0.8 

1.73 

445 

0.1 

60  ! 

Epsom  salts 

M.  S04-  7H,0 

0.8 

1.64 

445 

8 

0. 1 

50 

Washing  soda 

Na,  CO,*  10  H,0 

1.0 

1.44 

40C 

Table  salt 

N  Cl 

1.0 

2. 17 

285 

• 

0.44 

120 

Mahogany 

1.0 

255 

0.  32 

53 

Zonolite 

1.0 

200 

Cork 

1.0 

98 

0.6 

44 

0.44 

41 

0.2 

30 

Magnesium 

carbonate 

M^CO,*5H,0 

0.2 

1.69-.73 

68 

*  Taken  from  Hodgman  (1952). 
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The  first  attempts  to  use  the  camera  produced  reliable  values  only  for  the  larger 
(0.  5  to  1.  0  mm)  particles.  The  0. 1  mm  particles  were  not  identifiable  on  the  film. 
However,  this  method  was  used  on  a  number  of  test  materials  and  the  results  were 
extrapolated  in  an  attempt  to  narrow  the  investigation  to  a  few  materials  (Table  IV, 

Fig.  15).  By  means  of  an  auxiliary  closeup  lens  and  a  refined  focusing  technique,  it 
later  became  possible  to  identify  the  particles  in  the  0.  1  mm  sixe.  The  results  of 
the  free-fall  velocity  tests  are  shown  in  Figure  16  and  Table  IV. 

Coetficient  of  restitution;  Saltation  is  also  one  of  the  important  factors  in  the 
occurrence  of  blowing  snow  as  sKown  by  the  theoretical  study  presented  above.  It  is 
necessary  to  consider  in  detail  the  rebounding  of  snow  grains  when  they  strike  the  snow 
surface.  A  snow  grain  which  strikes  the  surface  during  the  saltation  process  cannot 
rebound  with  more  kinetic  energy  than  when  it  struck  (barring  outside  influences).  In 
practice  most  materials  will  rebound  with  some  lesser  amount  of  kinetic  energy.  Since 
the  mass  of  the  particle  remains  unchanged,  this  loss  of  kinetic  energy  can  be  expressed 
by  the  ratio  of  the  rebound  velocity  to  the  striking  velocity.  This  ratio  is  known  as  the 
coefficient  of  restitution  (see  p.  %  ).  According  to  Eshbach  (1936)  this  factor  depends 
only  on  the  materials  involved  in  the  impact  and  in  not  affected  by  their  size  or  shape, 
etc.  The  coefficient  of  restitution  is  considered  to  be  an  important  factor  in  the  salta¬ 
tion  process.  For  this  reason  it  was  decided  to  attempt  to  measure  the  values  of  the 
coefficient  for  snow  grains  and  some  of  the  potential  model-snow  materials. 

The  coefficient  of  restitution  was  measured  by  use  of  the  drop  chamber  and  the 
high-speed  camera  to  photograph  small  pieces  of  the  test  material  dropping  on  to  a  flat 
surface  of  e  same  material.  The  falling  and  rebounding  velocities  of  die  material 
were  calc  .<ited  and  the  ratio  of  the  two  gave  the  coefficient  of  restitution.  These  tests 
were  made  on  four  materials,  two  of  which  have  known  values  for  the  coefficient: 

a.  Steel  -  3/32  in.  carbon  steel  bearing  balls  dropped  on  a  steel  plate. 

b. .  Glass  -  1.0  mm  glass  beads  dropped  on  1/4  in.  plate  glass. 

c.  Borax  -  small  crystals  of  borax  (approx.  2  mm)  dropped  on  a 
large  borax  crystal  oriented  so  that  a  flat  surface  was  normal  to 
the  drop  line. 

d.  Ice  grains  -  small  pieces  of  ice  (1.  0  mm)  prepared  by  freezing 
drops  of  water  and  dropped  on  a  block  of  ice. 

The  entire  experiment  on  ice  grains  was  conducted  in  a  refrigerated  room  in  which  the 
air  temperature  was  -21C.  The  results  of  these  tests  are  given  in  Table  V.  The  dis¬ 
crepancy  between  the  published  and  observed  values  for  steel  may  be  due  to  difference 
in  the  type  of  steel  used.  The  type  of  steel  used  to  obtain  the  published  value  is  not 
cited  in  the  literature. 

Threshold  velocity:  A  study  of  the  relationship  between  particle  size,  parti¬ 
cle  density,  and  threshold  velocity  was  undertaken  in  the  pilot  wind  tunnel,  which  has  a 
test  section  of  10  in.  by  20  in.  by  8  ft.  In  the  first  test  a  number  of  materials  of  differ¬ 
ent  densities  were  prepared  to  a  uniform  particle  size  of  approximately  0.  5  mm.  Each 
of  these  materials  was  exposed  in  the  pilot  tunnel  while  the  velocity  of  the  wind  was 
gradually  increased  to  the  threshold  value.  The  results  (Fig.  17)  demonstrate,  as 
would  be  expected,  that  the  threshold  velocity  increases  with  increasing  particle  densi¬ 
ty.  In  a  second  test,  uniform  samples  of  cork  were  prepared  in  approximately  the 
following  sizes:  0.1,  0.3,  0.5,  0.8,  and  1.  0  mm  in  diameter.  Again,  samples  of  each 
were  exposed  in  the  wind  tunnel  and  the  threshold  velocity  was  determined  for  each  size. 
The  results  (Fig.  18)  show  that  the  threshold  velocity  diminished  with  decrease  in 
particle  diameter  down  to  a  minimum  at  approximately  0.  8  mm  diameter  and  then  in¬ 
creases.  The  relationship  between  threshold  velocity  and  particle  diameter  appears  to 
be  associated  with  the  viscous  surface  layer  of  air  into  which  very  small  particles  sink 
and  are  undisturbed  by  moderate  winds.  However,  the  thickness  of  this  boundary  layer 
diminishes  with  increasing  velocity  and  eventually  exposes  small  protuberances  and 
ripples  to  the  effect  of  higher  winds. 


33 


SCALE  MODEL  STUDIES  ON  SNOW  DRIFTING 


Selection  of  a  mow  aimu- 
Utor  from  the  meterUIe  teeted; 
From  the  reeulte  theee  etudiee 
it  was  decided  to  uee  commercial 
borax  (NaiB^Oy  ■  1  OH^O)  ae  the  snow 
simulator  in  die  first  series  of  tests 
in  the  large  wind  tunnel.  The  borax 
appeared  to  more  nearly  approach 
the  requirements  for  a  modeled  enow 
than  any  of  the  other  materials 
tested;  In  addition,  it  was  commer* 
cially  available  in  large  quantities, 
economical,  and  easily  handled. 

The  particle  else  distribution  of  a 
sample  of  the  borax  used  in  the 
wind  tunnel  is  shown  in  Table  VI. 

The  mean  sixe  of  the  particles  was 
larger  than  desired,  but  the  cost  of 
obtaining  it  in  0. 1  mm  size  was  pro¬ 
hibitive.  The  distribution  charac¬ 
teristics,  however,  compare 
favorably  with  those  for  snow  grains 
given  by  Diunin  (1956). 

A  photomicrograph  of  the  borax 
crystals  (Fig.  19)  shows  them  to  be 
quite  similar  to  old  snow  grains. 

WIND  TUNNEL  MODIFICATION 

Two  of  the  New  York  University 
wind  tunnels  were  used  in  various 
stages  of  this  project.  The  10  x  20 
in.  wind  tunnel  was  used  to  study 
snow  simulating  materials,  make 
preliminary  snow-drifting  experi¬ 
ments,  and  test  proposed  design 
changes  for  the  3^  x  7  ft  wind  tunnel 
where  the  final  snow-drifting  ex¬ 
periments  were  made. 

The  10  X  20  in.  wind  tunnel  has 
a  test  section  10  in.  wide,  20  in. 
high,  and  8  ft  long  in  the  direction 
of  air  movement.  One  wall  and  the 
ceiling  are  made  of  glass.  The 
wind  tunnel  draws  air  from  and  ex¬ 
hausts  to  the  laboratory  atmosphere. 
The  particles  which  are  drawn  out 
of  the  test  section  are  collected  in  a 
cloth  bag  which  covers  the  exhaust. 

The  3^  X  7  ft  wind  tunnel  has  a 
test  section  3^  ft  high,  and  7  ft  wide 
and  a  usable  length  of  30  ft.  The 
general  layout  of  the  test  section  is 
shown  in  Figure  20.  The  air  supply 
is  taken  from  within  the  laboratory 
building  and  first  passes  through 
various  air  straightening  devices 
before  it  reached  the  test  section. 
The  air  is  drawn  through  the  wind 


Figure  16.  Relationship  between  particle 
diameter  and  free-fall  velocity  for  all  par¬ 
ticles  measured. 


Figure  17.  Relationship  between  parti¬ 
cle  density  and  fluid  threshold  velocity. 
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Figure  18.  Relationehip  between  par¬ 
ticle  diameter  and  fluid  threshold  ve¬ 
locity  for  cork. 


Figure  19.  Photomicrograph  of  0.  1  mm 
borax  crystals  (enlargement* 25x). 


Table  V.  Coefficient  of  restitution  measurements. 


Material  Published  value*  Exp  value 


Steel  0. 55  0. 66 

Gl»*«  0. 95  0. 95 

Borax  0. 55 

Ice  grains  0.  33 


*Data  taken  from  Eshbach  (1936) 


Table  VI.  Sise  distribution  of  borax  used  in  the  tunnel  tests. 


Sieve  no. 

Material  retained 

Sample  A 

Vol.  (cm*)  %  by  vol. 

Sample  B 

Weight  (g)  %  by  weight 

18  (1000  |i) 

0 

0 

0 

0 

20  (840  p) 

trace 

<0. 1 

0 

0 

35  (500  p) 

trace 

<0.1 

0 

0 

50  (297  p) 

130 

8.0 

59.4 

11.2 

140  (105  p) 

1415 

87.0 

450. 85 

84.9 

170  (88  p) 

40 

2.5 

18.70 

3.5 

greater  t^n  170 

40 

2.5 

2.15 

0.4 
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Figure  20.  Sketch  of  test  section  of  New  York  University  3^  x  7  ft  wind  tunnel. 

tunnel  by  a  variable  spee^l  exhaust  fan  located  downstream  of  the  test  section  and  is 
exhausted  through  the  roof. 

The  major  modification  made  to  the  wind  tunnel  for  the  snow-drifting  experiments 
was  the  installation  of  a  hopper  to  introduce  the  snow  material  through  the  celling  at 
the  upstream  end  of  the  test  section.  This  served  to  produce  the  etfect  of  falling  snow 
which  is  one  of  the  features  of  the  snow  environment  on  the  Greenland  Ice  Cap.  High 
velocity  alone  will  not  always  cause  blowing  snow  on  the  Ice  Cap  (e.g.  •  Hamilton,  1958). 
The  falling  snow  appears  to  loosen  particles  previously  deposited  on  the  surface  and  to 
initiate  saltation  at  a  lower  wind  speed  than  would  otherwise  be  required  to  start  move¬ 
ment. 

The  hopper  was  de'qigned  to  handle  powdered  borax  which  had  been  selected  for  the 
snow  simulator.  A  smul  hopper  was  first  tested  in  the  10  x  20  in.  wind  tunnel.  Clog¬ 
ging  and  bridging  of  the  borax  created  difficulties  which  required  modification  of  the 
original  design.  When  the  first  series  of  experiments  was  conducted  it  became  evident 
that  further  modifications  would'be  necessary  for  efficient  handling  of  the  several  tons 
of  borax  used  in  the  experiments.  During  the  course  of  an  experiment,  the  moving 
snow  material  collected  in  an  enlarged  portion  of  the  wind  tunnel  downstream  of  the 
test  section  ahead  of  the  exhaust  fan.  It  was  periodically  removed  and  transferred  to 
the  hopper. 
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SCALE  MODEL  EXPERIMENTS  ON  SNOW  DRIFTING 
CHARACTERISTICS  OF  VARIOUS  STRUCTURES 

The  final  ecale  model  experiments  were  conducted  in  the  3|  x  7  ft  wind  tunnel 
described  earlier.  No  attempt  whs  made  to  scale  the  velocity  profile  in  these  experi¬ 
ments  .  The  modeled  objects  were  placed  approximately  1 5  ft  downstream  of  the 
location  where  the  simulating  snow  was  dropped  into  the  airstream  from  the  hopper. 
Before  an  experiment  was  started,  the  floor  of  the  test  section  was  covered  wiA  a 
layer  of  simulating  snow  which  extended  about  13  ft  upstream  of  the  modeled  objects. 
Upstream  of  the  snow  there  was  7  ft  of  bare  test  section  floor  surface  which  joins  the 
wind  tunnel  contraction  cone. 

In  the  initial  tests  it  was  not  possible  to  make  any  exact  measurements  of  the  snow 
drifts  around  the  various  models.  However,  an  excellent  set  of  photographs  were  ob¬ 
tained  for  all  the  tests;  and  most  of  the  tests  were  also  recorded  in  time-lapfe  motion 
pictures.  The  time-lapse  pictures  were  taken  every  l|  sec,  so  that  a  period  of  1  hr 
and  40  min  could  be  recorded  on  a  100-ft  roll  of  film.  Inspection  of  the  photographs 
and  the  motion  pictures,  plus  visual  observation  during  the  tests,  produced  the  qualita¬ 
tive  results  presented  in  this  report.  All  of  the  tests  were  run  at  a  tunnel  velocity  of 
19  ft/ sec  (13  mph)  which  corresponds  to  an  atmospheric  velocity  of  41  mph.  At  this 
tunnel  speed  the  borax  particles  seemed  to  be  just  above  the  fluid  threshold  velocity. 

Only  a  few  borax  particles  were  in  motion  when  the  material  was  not  falling  from  the 
overhead  hopper,  while  the  whole  tunnel  was  filled  with  blowing  borax  when  the  hopper 
was  feeding  the  windstream.  During  this  latter  process  saltation  was  very  well  devel¬ 
oped  and  surface  creep  was  also  very  active,  as  shown  in  the  time-lapse  motion  picture. 
Five  different  model  situations  were  studied  and  are  described  below. 

Clements -panel  buildings  oriented  crosswind;  Four  models  of  Clements- 
panel  buildings,  frequently  used  in  Arctic  locations,  were  cons,t,ructed  to  represent  a 
prototype  structure  oO  ft  long,  25  ft  wide,  and  10  ft  high.  They  were  erected  on  pil¬ 
ings  which  represented  a  5  ft  elevation  from  the  surface  of  the  snow*.  This  first 
building  was  placed  about  15  ft  (scaled)  downwind  of  the  slot  wnere  the  snow,  was  enter¬ 
ing  the  windstream.  The  structures  were  oriented  with  their  long  axes  normal  to  the 
wind  and  with  the  distance  between  the  buildings  equivalent  to  a  scaled  25  ft.  Figures 
21  through  28  show  this  test  at  various  stages  during  drift  accumulation. 

The  snow  was  introduced  into  the  windstream  for  a  period  of  3  hr.  This  "snow¬ 
storm"  was  approximately  equal  to  all  the  blowing  snowstorms  occurring  during  1  yr 
in  certain  areas  on  the  Greenland  Ice  Cap.  After  this  phase,  the  test  was  continued 
for  another  2  hr  with  only  the  wind  blowing.  This  later  phase  represents  a  period  of 
high  wind  without  precipitation  in  Greenland.  Such  winds  cause  erosion  of  the  surface 
with  loss  of  snow  at  one  location  and  deposition  at  another.  This  phase  of  the  tests  was 
undertaken  to  determine  the  interrelationship  between  erosion  and  building  or  structure 
design.  During  the  first  3  hr  the  structures  became  almost  completely  drifted  over 
with  snow.  However,  during  the  later  part  of  the  test,  a  considerable  amount  of  the 
snow  was  removed  from  around  the  structures  by  erosion.  The  photographs  give  a 
good  qualitative  illustration  of  the  processes  during  the  accumulation  and  erosion 
periods. 

Outstanding  in  Figure  21  is  the  stream  of  snow  blowing  off  the  roof  of  the  upwind 
structure.  The  stream  is  striking  the  side  of  the  next  building  and  most  of  the  snow  is 
falling  into  the  space  between  the  buildings.  With  this  layout  the  open  space  between 
buildings  could  be  expected  to  become  drifted  in  with  snow  very  rapidly.  In  Figure  22 
the  rocd  of  building  1  is  well  covered  with  snow  while  the  roofs  of  the  following  structures 
are  entirely  clear.  This  is  most  likely  due  to  a  combination  of  two  factors.  The  alter¬ 
ation  in  the  wind  velocity  pattern  due  to  the  blocking  effect  of  the  first  structure  and  the 
consequent  change  in  the  capacity  of  the  air  to  support  the  mass  of  snow  causes  signi¬ 
ficant  amounts  ctf  snow  to  be  deposited  in  front  of  and  on  the  roof  of  the  first  structure. 
The  second  structure  is  in  the  wake  of  the  first  and  the  air  passing  immediately  over  it 
is  unsaturated  with  respect  to  its  capacity  to  hold  snow  so  that  very  little  snow  accumu- 
lates  on  the  building  and  in  the  space  between  buildings  1  and  2.  This  is  substantiated 
*  'fhroughout  this  section  oi  this  report,  the  word  snow  is  used  to  signify  both  real  snow 
and  the  model  snow  (borax).  It  is  not  felt  that  any  confusion  will  arise  in  the  meaning 
whenever  it  is  used. 
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Figure  21.  Clement*  panel  building*  no*.  2  and  3 
(numbered  from  upwind)  oriented  cro**wind,  ^  hr 
after  atart  of  teat.  Wind  velocity  from  left  to  right. 


Figure  22.  Clement*  panel  building*  no*.  1,  2, 
and  3  oriented  croaawind,  1  hr  after  atart  of  teat. 
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Figure  23.  Side  view  of  buildings  nos.  1,  Z, 
and  3  at  the  same  time  as  shown  in  Figure  22. 


Figure  24.  Clements  panel  buildings  nos.  1,  2, 
and  3  oriented  crosswind,  3  hr  after  start  of  test. 


Figure  25.  Clements  panel  buildings  nos.  3  and 
4  oriented  crosswind,  3  hr  after  start  of  test. 


Figure  26.  Clements  panel  buildings  nos.  1,  Z,  and  3  oriented 
crosswind,  3|  hr  after  start  of  test,  |  hr  after  start  of  erosion. 
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Figure  27.  Schematic  view  of  drift 
pattern  after  3  hr  of  snow  deposition, 

by  Figure  23  which  clearly  shows  the  eddy  pattern  over  the  first  building.  The  free 
stream  pattern  lowers  after  passing  the  first  building  so  that  it  sweeps  the  roof  of  the 
next  building  clear. 

Figures  24  and  25  show  the  buildings  at  the  end  of  the  snow-storm  period.  It  can 
be  seen(Fig.  27)  that  the  final  accumulation  pattern  is  related  to  the  events  which  oc¬ 
curred  during  the  testing  period.  Building  1  is  completely  inundated  in  front  and  the 
roof  is  still  covered;  building  2  and  the  space  between  1  and  2  remain  relatively  clear. 
The  space  between  buildings  2  and  3  contains  a  drift  with  a  very  sharp  iqiwind  slope. 

This  agrees  with  the  flow  pattern  shown  in  Figures  21  and  23.  The  space  between 
buildings  3  and  4  has  been  completely  filled,  and  the  roofs  of  buildings  2,  3  and  4  have 
remained  clear.  The  results  of  this  phase  of  the  test  indicate  that  the  spacing  between 
buildings  may  have  been  insufficient  to  permit  high  velocity  flow  during  the  accumulation 
period.  The  build-up  of  snow  around  the  buildings  shown  in  Figure  26  appears  to  be  due 
in  part  to  the  coalescence  of  drifts  produced  by  each  building.  Roots  and  Swithinbank 
(1955)  reported  such  coalescence  of  drifts  at  Maudheim  in  the  Antarctic  where  spacing 
between  Duildings  was  inadequate. 

After  the  erosion  phase  began  (Fig.  26),  the  roof  of  building  1  started  to  clear  and 
the  space  between  buildings  1  Md  2  started  to  fill  up.  This  is  probably  explained  by  the 
fact  that  surface  creep  is  a  very  important  factor  in  this  phase.  After  2  hr  of  erosion 
(Fig.  28),  the  roof  of  building  1  is  cleared  and  the  entire  drift  pattern  is  decreasing  as 
the  strong  wind  constantly  removes  small  amounts  of  snow  from  the  test  area. 

Clements  pwel  builiUnas  oriented  downwind;  The  same  four  buildings  used  in 
the  first  test  were  also  tested  vnth  their  major  axes  parallel  to  the  windstream.  The 
separation  between  each  building  again  represented  a  scaled  distance  of  25  ft.  This 
arrangement  was  subjected  to  a  snow  deposition  period  of  3  hr  with  all  other  conditions 
maintained  similar  to  the  first  test.  Inspection  of  the  photographs  (Fig.  29,  30)  shows 
that  the  leading  edge  of  the  roofs  of  the  upwind  buildings  became  snow-covered,  while 
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in  the  central  area  the  roofs  were  essentially  clear.  This  was  a  consistent  feature  of 
the  tests  with  buildings  erected  with  the  long  axis  parallel  to  the  wind  flow.  The  ori¬ 
entation  of  the  buildings  with  their  long  axes  parallel  to  the  wind  Appeared  to  induce 
channelising  of  wind  flow  and  increase  velocity,  resulting  in  less  deposition  against 
the  side  walls  of  the  buildings  after  the  clear  under-space  was  saturated.  It  is  not 
known  whether  this  channelising  and  self-clearing  action  would  prevail  on  the  Green¬ 
land  Ice  Cap  where  winds,  although  predominantly  from  one  direction,  can  shift 
through  20  deg  during  the  progress  of  a  major  storm.  In  addition,  under  natural  con¬ 
ditions  storm  winds  accompanied  by  falling  and  drifting  snow  may  be  from  a  direction 
45  deg  divergent  from  clear  weather  winds.  Drifts  accumulated  between  buildings  by 
wind  from  one  direction  may  be  protected  from  erosion  by  die  buildings  when  the  wind 
shifts  a  few  degrees. 

Elevated  radar  type  building;  A  model  was  constructed  (Fig.  31,  32)  to  repre¬ 
sent  a  building  sO  ft  wide,  60  ft  long7  and  10  ft  high  set  on  columns  which  gave  it  a  10-4t 
elevation  above  the  snow  surface.  This  model  was  exposed  to  3  hr  of  snow  deposition. 
The  most  significant  feature  of  this  test  was  that  the  blowing  snowstorm  initially  caused 
erosion  under  the  model.  This  erosion  was  great  enough  to  remove  the  entire  layer  of 
snow  under  the  model  and  reveal  the  floor  beneath.  This  test  was  repeated  with  a 
much  deeper  layer  of  snow  under  the  model.  Again  erosion  occurred  to  such  an  extent 
that  the  stability  of  the  model  was  threatened. 

This  same  model  was  modified  so  that  it  was  supported  on  a  buried  foundation  or 
crib  (Fig.  33).  The  model  foundation  used  was  a  rectangular  frame  of  2  by  4  in.  beams 
buried  flush  ^th  the  snow  surface.  Erosion  began  at  the  start  of  the  test  but  quickly 
reached  an  equilibrium  pattern  shown  by  the  dashed  line  (Fig.  33).  At  no  time  was  the 
stability  of  the  structure  threatened.  After  equilibrium  was  reached,  drifting  proceeded 
in  the  usual  manner  (Fig.  34).  As  a  corollary  to  this  experiment,  a  model  was  buried 
so  that  its  flat  roof  was  flush  with  the  snow  surface.  No  drift  snow  accumulated  on  the 
roof  or  in  the  vicinity  of  this  structure.  This  approach  to  the  problem  would  be  ideal 
except  that  many  structures  (particularly  radar  buildings)  must  have  above -surface  ex¬ 
posed  sections  to  house  electronic  detection  equipment. 

Snow  fence  mo^ls;.  Models  of  two  different  types  of  snow  fences  were  con¬ 
structed  an^  tested  briefly  to  see  if  the  3^  by  7  ft  tunnel  produced  a  pattern  which  was 
comparable  with  observed  natural  conditions.  The  first  of  these  was  a  10:1  scale  model 
of  one  of  the  drift  fences  used  by  the  Union  Pacific  Railroad.  This  fence  was  tested  in 
both  a  vertical  (Fig.  35)  and  a  60*  position  (Fig.  36).  The  other  fence  tested  was  a 
10:1  model  of  a  paper  strip  fence  (Fig.  37,  38)  previously  tested  on  the  Greenland  Ice 
Cap  and  discussed  by  Gerdel  (1960).  The  accumulation  of  snow  by  this  fence  during  the 
test  period  (3  hr)  agrees  in  a  qualitative  sense  with  the  pattern  stiown  by  Gerdel.  The 
accumulation  by  the  Union  Pacific  type  fence  also  agreed  with  tne  classical  atmospheric 
drift  pattern  for  both  angles  of  inclination.  Differentiation  between  the  characteristics 
of  the  two  different  angles  was  not  attempted  at  this  time. 

In  another  test  the  Union  Pacific  type  fence  was  set  up  at  a  90*  angle  with  the  V  of 
the  angle  pointing  upwind  (Fig.  39).  This  arrangement  produced  a  very  long  area  of 
minimum  accumulation  downwind  of  the  fence  (Fig.  40).  The  length  of  this  area  was 
equal  to  approximately  25  times  the  height  of  the  fence.  Since  the  fence  represented 
an  actual  height  of  8  ft,  this  clear  area  represents  an  extent  of  about  200  ft.  This 
fence  shape  might  prove  very  effective  in  keeping  a  moderate -siaed  area  free  from 
drifting.  Further  testing  of  this  arrangement  appears  to  be  very  desirable. 

Fuel  storage  Mn^  In  die  last  test  a  number  of  cylindrical  tanks  were  pre¬ 
pared  which  represented  POL  storage  tanks  with  a  hsight  od  10  ft  and  a  diameter  of 
ft  (Fig.  41,  42).  These  were  subjected  to  snow  deposition  in  the  same  way  as  the 
oUier  models.  These  tanks  developed  the  same  trough  pattern  on  ttieir  downwind  side 
as  did  the  V-shaped  fence.  The  base  of  each  tank  remained  clear  for  a  period  of  time, 
but  when  the  snow  near  the  tanks  became  high  enough,  it  began  to  slide  down  into  the 
empty  space,  eventually  filling  it  up. 
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Figure  28.  Clemente  panel  buildings  nos.  1,  2,  and  3  oriented 
crosewind,  5  hr  after  start  of  tests,  2  hr  after  start  of  erosion. 

The  borax  has  a  lesser  angle  of  repose  and  less  coherence  than  snowt  consequently 
it  may  be  presumed  that  the  vortex-produced  clearing  around  each  tank  may  have  been 
less  deep  in  the  wind  tunnel  test  dian  it  would  be  under  natural  conditions  unless  the 
natural  snow  is  sufficiently  hard  to  resist  erosion. 

FUTURE  WORK 

Further  investigation  in  this  field  of  synthesised  blowing  snow  environments  is 
needed.  The  theoretical  study  of  the  blovring  snow  process  should  be  given  more  empha¬ 
sis.  The  correlation  between  the  theoretical  work  and  the  actual  movement  of  synthetic 
snow  particles  should  be  investigated  in  more  detail. 

Much  further  work  also  can  be  done  on  the  10:1  scale  testing.  All  tests  should  be 
repeated  with  emphasis  on  the  quantitative  aspects  of  the  drifts.  In  addition  many  differ¬ 
ent  types  of  snow  fences  and  other  structures  should  be  tested. 

A  program  should  be  initiated  to  develop  a  synthetic  snow  material  which  could  be 
used  for  a  much  smaller  scale,  e.g.  50:1.  This  would  greatly  increase  the  number  and 
sine  of  models  which  could  be  studied  in  the  wind  tunnel  and  permit  investigation  on  the 
interaction  between  units  in  large  installations  conqposed  of  many  buildings  and  different 
types  of  structures. 
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Flgur*  29.  Clement*  panel  building*  oriented  with 
the  wind,  2  hr  after  *tart  of  ta*t,  looking  downwind. 


Figure  SO.  Clemonta  panel  building*  oriented 
with  the  wind,  2  hr  after  start  of  teat.  Sid*  view. 
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Figure  31.  Radar  type  building,  1  hr  after  atart  of  teat. 


Figure  32.  Radar  type  building,  1^  hr  after  atert  of  teat. 
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Figure  33.  Schematic  view  of  foundation  used 
to  support  radar  type  building  for  test  rerun. 


Figure  34.  Radar  type  building  with  spe¬ 
cial  foundation,  2  hr  after  start  of  test. 
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Figure  35.  Union  Pecific  type  enow  fence 
inclined  90*,  3  hr  after  start  of  test. 


Figure  36.  Union  Pacific  type  snomr  fence 
inclined  60*,  Z  hr  after  start  of  test. 
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Fifur*  37.  Modal  of  papar  atrip  fanea 
tastad  on  tha  Graanlatad  lea  Cap. 


Fifttra  3S.  Papar  atrip  fanea. 
3  hr  aftar  atart  of  taat. 


Figure  39.  V-shaped  Union  Pacific 
type  fence  prior  to  start  of  test. 
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Figure  40.  Clear  area  downwind  of  V-shaped 
fence,  after  2-hr  test  period. 
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Figure  42.  Fuel  tanks  after  3  hr 
of  blowing  snow,  looking  downwind. 
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